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ABSTRACT 



The following describes a model of a broad planetary ring whose sharp edge 
is confined by a satellit e's rrfi^ Lindblad resonance (LR). This model uses the 
streamline formalism of iBorderies et al.l (119821 . Il985l ) to calculate the ring's in- 
ternal forces, namely, ring gravity, pressure, and viscosity. The model also allows 
for the possibility of a drag force that can affect small ring particles directly, and 
large ring particles indirectly via collisions with the small. The model calculates 
the streamlines' forced eccentricities e, their longitudes of peripase cD, and the 
surface density a throughout the perturbed ring. This model is then applied 
to the outer edge of Saturn's B ring, which is maintained by an m = 2 inner 
LR with the satellite Mimas. A suite of ring models are used to illustrate how 
a ring's perturbed state depends on the ring's physical properties: its surface 
density, its viscosity, the ring particles' dispersion velocity, and the strength of 
the hypothetical drag force. A comparison of model results to the outer B ring's 
observed properties suggests that the ring's surface density there is 10 < a < 280 
gm/cm^ in the ring's outermost ~ 40 km. 

The ring's sharp edge identifies the site where the ring's viscous torque pre- 
cisely counterbalances the perturbing satellite's gravitational torque on the ring. 
However, an examination of several seemingly conventional viscous B ring models 
shows that they all fail, by wide margins, to balance these torques at the ring's 
outer edge. This is partly due to the ring's self-gravity, which tends to reduce 
forced eccentricities near the resonance. But this is also due to the fact that a 
viscous ring tend to be nearly peri-aligned with the satellite. Both effects con- 
spire to reduce the satellite's torque on the ring, which in turn makes the ring's 
edge more difficult to maintain. Nonetheless, the following shows that a torque 
balance can still be achieved in a viscous B ring, but only in an extreme case 
where the ratio of the ring's bulk/shear viscosities satisfy Vb/i^s ~ 10^. However, 
if the dissipation of the ring's forced motions is instead dominated by a weak 
drag force, then the satellite can exert a much stronger torque across a wider an- 
nulus in the ring, which can successfully counterbalance the ring's viscous torque 
there. We also show how this streamline model can be adapted to study other 
interesting ring phenomena, such as narrow eccentric ringlets and nonlinear spiral 
density waves. 



Subject headings: planets: rings 
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1. Introduction 

The outer edge of Saturn's main B ring is confined by an m = 2 inner Lindblad resonance 
(ILR) with the satelhte Mimas, while the outer edge of the main A ring is confined by m = 7 
ILRs with the coorbital satelhtes Janus and Epimetheus (Porco et al 1984). Ring particles 
orbiting near a resonance execute m radial excursions, or epicycles, which gives a ring's edge a 
scalloped, m-lobed appearance. The ring particles' epicyclic amplitude Rm, which is the half- 
amplitude of the particles' radial excursions, is obviously governed in part by the mass of the 
perturbing satellite. However the ring's internal forces — self gravity, pressure, and viscosity — 
also play a role here by tempering the ring's response to the satellite's resonant gravitational 
perturbations. Consequently, modeling these rings in a way that accounts for all of these 
perturbations, and then comparing those models to observations of the ring's edge, should 
allow one to assess the relative importance of a ring's various internal forces. This activity will 
also allow one to make estimates of, or place limits on, the ring's physical parameters, such 
as its surface density a, viscosity u, and the ring particles' dispersion velocity c. Also note 
that small ring particles are susceptible to drag forces — plasma drag, Poynting-Robertson 
drag, and/or the Yarkovsky effect — whose significance can also be assessed by this kind of 
modeling, too. 

The epicyclic amplitudes of the outer A and B rings are quite small, ~ lO's of 
km, while the ring's radii are r ~ 10^ km, so the ring particles' noncircular motions are 
only ~ 0.01% in a fractional sense. The smallness of those motions also makes any time- 
dependent ring models, such as Nbody, SPH, or hydrodynamic treatments, very difficult, 
due to the very many particles or hydrodynamic cells needed to simulate the ring-edge's 
entire circumference. Also keep in mind that one does not know the ring's equilibrium 
state in advance, so simulations of a ring's time-evolution would initially be dominated by 
transient motions that can easily mask the ring's small forced motions. Consequently, time- 
dependent models must also evolve the system until the ring's dissipation has damped out 
those transients, which can also cost a lot of CPU time. 

Semi-analytic methods instead provide a much more efficient and illuminating method 
for studying perturbed planetary rings. These methods are enabled by the streamline con- 
cept, which is a very powerful tool for studying planetary rings. A streamline is the epicyclic 
path that is traced by numerous ring particles that all have a common semimajor axis. A 
planetary ring can then be thought of as a sum of many such streamlines. This streamline 
approach also allows one to calculate the forces that any one streamline exerts on any one 
ring particle. Summing over all the forces exerted by all of the ring's streamlines, and in- 
serting those forces into Newton's equations of motion, then provides a model that can be 
used to calculate the perturbed ring's equilibrium state. 
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Streamlines were first used to show how a narrow ringlet's self- gravity can counter the 
differ ential precession that occurs when orbiting an oblate planet (iGoldreich &: Tremaine 
19791 ). A streamline approa ch was also used to sh ow how viscosity helps a satellite maintain a 
planetary ring's sharp edge (iBorderies et al.lll982l ). and streamlines were used to examine how 



the g ap-embedded satellite Pan maintains the Encke gap in Saturn's A ring ( Borderies et al. 



1989 ) . Streamlines are also used to study nonlinear spiral density waves (jBorderies et al. 
19861 : iLongaretti fc Borderieslll986l ). Evidently, the streamline concept is a very general tool 
that can be used to study a variety of interesting ring phenomena. And in the following, we 
use a streamline formalism to examine in detail how a satellite's m**^ ILR can disturb as well 
as maintain the sharp outer edges of Saturn's main rings. 

Section [2] reviews the subject in detail, and derives all of the equations that will be 
needed to assemble a comprehensive model of a broad planetary ring whose edge is confined 
by a satellite's m^^ Lindblad resonance. However, a reader who is uninterested in the many 
details and derivations can skip ahead to Section [3l which examines a suite of B ring models 
that illustrate how the perturbed ring's state — its epicyclic amphtude, its orientation, and 
it surface density variations — all depend on the ring's physical properties: a, u, and c. This 
section also shows how observations can be used to infer a ring's physical properties, which 
for the B ring are quite unknown. Section H] then discusses some important side issues, while 
Section [5] summarizes our main findings. 



2. Equation of motion 

This Section derives the equation of motions that will be used to calculate the motion 
of an orbiting ring particle while it is perturbed by an orbiting satellite and the planetary 
ring. The particle's position vector r{t) evolves over time t according to Newton's second 
law of motion, 

where $(r) = $p + $s is the total gravitational potential that is due to the central planet $p 
and satelhte and a is the acceleration of the particle that is due to the forces exerted by the 
planetary ring. For simplicity the following assumes that the satellite's orbit is circular and 
coplanar with the ring plane, but these results are unchanged if the satellite's orbit is instead 
slightly noncircular or inclined. Polar coordinates will also be used, where r = (r, 6) with 
r{t) being the particle's distance from the planet's center, and 6{t) its longitude measured 
relative to some x axis, with the z axis perpendicular to the orbital plane. Solutions to the 
equation of motion ([T]) are then obtained after Fourier expanding the perturbations that are 
acting on the orbiting ring particle. 
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2.1. Fourier expansions of the perturbations 

Planet-centered coordinates will be used, so the satellite's gravitational potential has 
direct and indirect parts that are 

r, 9) = + -^r ■ r„ 2 

|r — rs| 

where G is the gravitational constant, is the satellite's mass, and = {as,6s) is its 
position vector in polar coordinates, with being the satellite's semimajor axis and 6s its 
longitude. A Fourier expansion of that potential is 

^s{r,e) = ^0°(r) + 5Je ( C(Oe^'"^'-'=M (3) 



\m=l 



where the potential's Fourier amplitudes 0™ are 

<PTir) = - r $s(r, ^) cos{m^)d^ = - (36mi] (4) 

where cp = 6 — 6s is the particle's longitude relative to the satellite's, and the Laplace 
coefficient 

n Jq (1 + /^^ - 2/^ cos 

is a function of the ratio (3 = r /as. The Kronecker delta 5mi in Eqn. (jlj) is due to the indirect 
part of the potential in Eqn. ([2]), which only contributes to the m = 1 part of the satellite's 
gravity. 

The acceleration that the ring exerts on the particle is a(r, 6) = arV + agd, where ar and 
ae are the radial and tangential components. A Fourier expansion of those accelerations will 
also have the form 

ar{r,9) = A0(r) + 3fte ||^^A™(r)e*'"(^-^»)j (6a) 
af){r,9) = A°(r) + sRe [ ^A^(r)e''"(^-^=M . (6b) 



These Fourier expansions are convenient since each of the m > 1 terms Eqn. ([3]) correspond 
to Lindblad resonances that are all spatially segregated. Consequently, when solving the 
equation of motion for the particle's motion, we only need to retain a single m^^ term in the 
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expression for $s. This is also true for the ring's internal accelerations and ag, since they 
are excited by the satellite's m^^ resonant perturbation of the ring. In light of this, write 

$,(r,e) ^ 0^(r)e*'"(^-^^) (7a) 
ar{r,e) ^ A°(r) + A;"(r)e^"(''-''^) (7b) 
ae{r,e) ~ A°(r) + A^(r)e^"^(''-^^) (7c) 

with the 3fJe() notation dropped henceforth, so it is to be understood that one is to preserve 
only the real parts of the following equations. Also note that the axisymmetric part of the 
satellite's potential, 0°, was omitted from Eqn. (I7aj) . since it is convenient to combine it with 
the planet's potential, $p ^ $p + (j)^, which is also axisymmetric. Lastly, note that 0™, A^, 
and Ag are all real, while the AJT* and A"^ can be complex. 



2.2. motion near a Lindblad resonance 

The radial and tangential parts of the particle's equation of motion are 

r - rr = — 7^ + ar (8a) 
or 

r at r o9 



when Eqn. ([7]) is inserted into Eqn. ([T]). Similar equations are solved in Goldreich &: Tremainel 



( 119821 ) for an isolated particle that does not experience a perturbation from the ring {i.e. 



a = 0), so the following solution to the more general sl ^ pro blem will use a strategy and 



notation similar to that given in I Goldreich &: Tremainel (119821 ). A ring particle's orbit will 



be nearly circular, so its trajectory has the form 

r(t) =ro + ri(t) and 6{t) = 9o + Qot + 9i{t) (9) 

where the constant tq is the particle's mean distance from the planet, Qq is its mean angular 
velocity about the planet, and is an arbitrary phase. Since the perturbing accelerations 
a and — V$s are all small compared to the central planet's gravity — V$p, the particle's 
displacements from a purely circular orbit will be small such that |ri| <^ ro and <^ 1, 
which then allows the equation of motion to be linearized. 

The satellite's orbital angular velocity is ^2^, so its longitude is Og = Qgt when time 
if: = is chosen to be the time when it traverses the x axis. This also means that the 
particle's relative longitude that appears in Eqn. ([7]) is m{6 — 6s) = m9o +uJmt + mO\ where 
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i^m.{r) = m{Qo — Qs) is the particle's Doppler-shifted forcing frequency. But this quantity is 
usually needed only to lowest order in the small angles, so 

m{e - Os) ~ mOo + Umt. (10) 

The particle's specific angular momentum is /i = r^6', and the time-evolution of that 
quantity is obtained from Eqn. (]8bp with 9 = h/r"^: 

^ = _^ + rae = -zmCe''"^'"'^^ + rA^e + rA^e'^^'-'^\ (11) 

The first term is the specific torque that the satellite exerts on the ring particle, while the 
other terms are the specific torques that are due to the ring's internal forces. Evidently, 
the total torque dh/dt is the sum of secular {i.e., non-oscillatory) terms like tqA^q plus other 
oscillatory terms. In light of this, write h{t) = hQ{t) + hi{t) where = TqQq is the secular 
part of the particle's specific angular momentum h{t), while hi{t) is the oscillatory part. 
The secular and oscillatory torques on the particle, dho/dt and dhi/dt, are then 

"^^^ ~ roA° + T. (12a) 



dt 
dhi 



[-imcl)'^ + rov4^)e'(^"+"™*) (12b) 



when Eq. (11 II) is written to lowest order in the particle's coordinates, i.e.., with r ^ tq 
and m{9 — Og) — mOo + Umt. Note that we have also added to Eqn. f ll2ap an additional 
term Tg to represent the secular part of the specific torque that the satellite exerts on the 
particle; Section 12.61 will show that this second-order term is important only near the ring's 
sharp edge. And if the ring particle's orbit is to be static such that its mean orbit radius 
ro is constant, then the secular torque on the particle, dho/dt, must be zero, for otherwise 
that particle (as well as its neighboring ring particles) would drift radially. Consequently, 
equilibrium thus requires all particles to satisfy the torque-balance equation, Tg = —roAg. 



Integrating Eqn. fll2bp with respect to time t provides the oscillatory part of the part 



cle's specific angular momentum. 



Also note that \hi\ is small compared to = t^Qq. The quantity hi is then used to solve 
the radial part of the equation of motion, which is 

r-^ = + A^.+ (-^ + aA e^irnOM) ^^^^ 
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when 9 = h/r"^ and Eqns. ([7]) are substituted into Eqn. fl8a|) . Inserting r = tq + ri and 
h = Hq + hi into the above, Taylor-expanding to first order in the small quantities ri and 
hi, and then inserting h^ = t^Qq and Eqn. f[T^ into that result then yields 

( t - - (^^^^ - ^) - - 

(15) 

where it is understood that all quantities in the above are to be evaluated at r = tq. 

The ring particle is assumed to be in torque balance, so tq is constant and tq = 0. And 
since the terms in the first set of parentheses in Eqn. f[T^ are secular while the remaining 
terms are oscillatory, that parentheses, which is the condition for centrifugal equilibrium, 
must separately sum to zero, which provides the particle's mean angular velocity Qq = fl{ro) 
where 

= - ^. (16) 

The constant in the second parentheses in Eqn. ( fT5|) is the particle's epicyclic frequency, 
kq = ^^('"o); where 

-'-sn'^^. (17) 

Since the ring's radial acceleration \A^\ is small compared to the central planet's gravity, 
these angular frequencies are 



n ^ n,[i-^j (18a) 

and ^ « 11-^] (18b) 



where f2p = ^/7^~^'d^^^Jdr would be the particle's angular velocity if ring forces were absent, 
and Kp = ^ AVtp + rdVL^ / dr would be its epicyclic frequency when = 0. 

The terms on the right of Eqn. f[T^ that involve 0™ are the satellite's forcing function, 

or rUJrn 

which accounts for the satellite's radial and tangential forcings. All the coefficients on the 
right hand side of Eqn. (fT5l) will be known as the system's complex forcing function, 

■qf^{r) = ^T- + (20) 

Inserting this into the above then casts Eqn. (fT5|) in its simplest form, 

fi + Kin ~ ^^(ro)e*("^^o+^™*). (21) 



- 9 - 



2.3. single particle motion 

When the ring's internal forces are absent, i.e., a = 0, then \E'™ = is a constant, 
and Eqn. fl2T|) describes a driven simple harmonic oscillator whose solution is 

ri(t) = -i?„e^(™^»+^™*). (22) 



This solution is examined in iGoldreich fc Tremaind (119821 ). which is summarized here since 



those results are used throughout this study. Inserting Eqn. (1221) into Eqn. (1211) then yields 
the particle's epicyclic amplitude i?^, which is 

= (23a) 

where D{r) = - cj^ (23b) 

is the particle's distance from resonance in frequency- squared units. When a particle is far 
from a resonance, \D\ is of order Q^, and the particle's epicyclic amplitude Rm is negligibly 
small. However, when the particle is near a resonance, \D\ ^ and the particle's response 
to the satellite's resonant forcing is much larger. Exact resonance is the site where D{r^) = 0, 
or where k = ecj^ with e ± 1. If the central planet's potential $p were Keplerian, then 
$p = —GMp/r where Mp is the planet's mass and k = n = ^yGMp/r^. Inserting this 
into the resonance condition k = eUm = em{Q — Qg) then yields the radius of the Lindblad 
resonance, 

rr = (l- tjmfl^as (24) 

where is the satellite's orbit radius. Resonances having e = +1 are inner Lindblad 
resonances (ILRs) since they reside interior to the satellite's orbit, while those with e = — 1 
are outer Lindblad resonances (OLRs). The focus of this work will be on the outer edge 
of Saturn's main B ring, which is confined by an ILR with a satellite that orbits exterior 
to the ring, so e = +1 here. Lastly, note that if the central planet is oblate, then the 
particle's epicyclic frequency k differs slightly from its angular velocity fi, so Eqn. (12^ is 
only approximately true. However, a precise calculation of the resonance location is also 
given in Section 12.3. 1[ 

For a particle orbiting near a LR it is sufficient to linearize D(t) via D ^ xD where 
X = [r — rj)/rr is the particle's fractional distance from resonance, and 

dD 



V = r- 

dr 



3e(m - e)nl. (25) 



Inserting Eqn. (jl]) into (1T9|) and noting that ujrn — near a LR also provides the satellite's 
forcing function, which is 

^r(r.) = efTlisTo^l (26) 
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where Hs is the satelhte's mass in units of the central planet's, and /™ is 



^ 8(3 



(27) 



which is a positive numerical coefficient that depends on the resonance in question. For 
instance, this study is interested in Mimas' m = 2 ILR in the B ring, which has e = +1 and 
/J" = 1.500. Inserting these quantities into Eqn. fl23ap then provides the particle's epicyclic 
amplitude. 



3(m — e)x' 



as well as its forced eccentricity 



where tps = 



\ri\ _ 
ro 



3(m — e)\x\ 



'4's 



X 



\roV\ 3{m — e) 
is a dimensionless version of the satellite's forcing function. 



(28) 

(29a) 
(29b) 



Note that Rm is real and has the opposite sign as x. This then allows a parametrization 
of the particle's radial displacement ri in terms of its longitude 9, 



ri{9) ~ sgn(x)|i?m| cosm(6' — Og 



(30) 



when Eqn. (fTOj) is utilized. Also note that the magnitude of the particle's radial excursion 
|ri| is maximal when it arrives at that satellite's longitude, 6 = 9s- Consequently, if the 
particle's orbit lies interior to the resonance with x < 0, then ri < 0, so the particle must 
also be at periapse since r{6 = Og) = tq + ri = ro — \Rm\- Similarly, a particle orbiting 
exterior to resonance would be at apoapse, r = ro + \Rm\i when in conjunction with the 
satellite. So in summary, orbits interior to a LR are expected to be peri-aligned with the 
satellite's longitude, while orbits exterior should be apo-aligned. 

Interestingly, the outer edge of Saturn's B ring does not behave according to these 
expectations. Voyager observati ons revealed that the B ring's outer edge actually lies about 
24km exterior to the resonance (jPorco et al.lll984l ). This is quite a curiosity, since, if Mimas' 
m = 2 ILR is indeed responsible for maintaining this ring's outer edge, then one might expect 
the ring's outer edge to lie at or just interior to the resonance. Another curiosity is that 
the ring is peri-aligned with Mimas, where as one woul d expe ct the ring to be apo-aligned 
if it truly extended beyond the resonance. iPorco et al.l (Il984 ) suggest that this interesting 
behavior might be a consequence of the ring's internal forces, which are considered in Section 
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2.3.1. resonance location 



A Lindblad resonance is the site r = that satisfies D{rr) = 0, which is also where 

K,{rr) = eujm{f^r) = em[Q{rr) — ^Ig], (31) 

with the frequencies Q and k defined by Eqns. f|T6l) and f|T7|) . Those quantities also depend 
on gradients in the central planet's gravitational potential $„, which is 



$p(r) 



GM„ 



k=i \ py 



-2k 



(32) 



for bodies orbiting in an oblate planet's equatorial plan e, where the J2k are the planet's 
zonal harmonics and the are Legendre polynomials (IMurray &: DermottI Il999l ). All of 



the Saturnian system's physical constants that are used in the model, GMp, J2k, etc., are 
listed in Table [H For instance, solving Eqn. ( 13T1) for the location of Mimas' m = 2 ILR yields 
Qr = 117, 553.71 km when the standard Rp = 60, 330 km is adopted as Saturn's radius. The 
position Oj. will be referred to as the nominal resonance position, since exact resonance can 
be shifted further by the ring's internal forces, which are assessed below. 



2.4. the ring's internal forces 



A ring particle is also subject to additional forces, such as the gravity exerted by the 
entire ring, and collisions with other ring particles. The cu mulative effect of those collis ions 
will be modeled here as if they were due to pressure (c.f., iGoldreich fc Tremainelll978l ). A 
ring parti cle is also subiect to viscous effects, which could be due to collisions among ring 
particles (IGoldreich fc Tremaind Il982l ). or perhaps due to transie nt self-gravitatin g wakes 
in the ring (jPaisaka et al.ll200ll ) and/or the viscous overstability (ISalo et al.l 120011). Small 
ring particles are also susceptible to a number of other drag forces (IGoldreich fc Tremaine 
19821 ) that can be be accounted for here, too. The following subsections will assess the 
ring's internal forces, t he ca lculation of which is greatly facilitated by the streamline con- 
cept. iBorderies et al.l (119821 ) note that Eqn. (I5U]) also represents a streamline in the ring, 
which is a closed trajectory that is traced by numerous ring particles having a common 
semimajor axis tq. A broader planetary ring can then be re garded as being compr ised of 
numerous streamlines having distinct semimajor axes tq. As IBorderies et al.l (119851 ) show, 
the calculation of the acceleration a that the entire ring exerts on a single particle can be 
simplified substantially by considering the particle's interaction with individual streamlines. 
But doing that will first require understanding the relationship between the ring's surface 
density a{r, 6) and a streamline's semimajor axis a = tq, eccentricity e, and its longitude of 
periapse u. 
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2.4- 1- streamlines 

Equation fl20l) shows that if the ring's internal forces have a tangential component, then 
the equation of motion fl2Tl) for a ring particle's radial displacement ri is complex, so write 

ri(t) = _7^^e^(™^«+^'"*) ~ -7^^e^"^(''~''^) = -i^^e^^^^"^^"'^) (33) 

where 7^m(a) = -Rm.e"'""'^ is the particle's (or streamline's) complex epicyclic amplitude, 
and uj{a) its longitude of periapse, both of which are to be regarded as functions of the 
streamline's semimajor axis tq = a. The particle's total planetocentric distance is then 

r{a,e) = a + ri ^a-3^e[7^„(a)e^'^(^-''^)] = a - Rm{a) cosm{e - 6,, - Co). (34) 

Also identify e(a) = Rm/o. as the particle's forced eccentricity, with e and Rm to be regarded 
as positive real quantities. 

The following will also need the ring particle's longitudinal displacement, 6i, which is 
obtained by linearizing its total specific angular momentum h = r'^6 ~ rgilo+'^o^i+S^^o^o^i = 
Hq + hi, so ho = VqQq as expected and 



. JJ}^ + '21 = ?^ _ ^ iill ) (35) 
which is integrated to obtain 

«-(-^^^^f)^'-- (-) 

since Um ~ eK ~ eil near a LR. However the particle's epicyclic amplitude TZm will be large 
near a resonance, dwarfing the other terms, so 

Ox ^ 37 

a 

and 

Qx ~ 2eJ^oexp \im[B - 6^ - Co)]. (38) 



2.4-2. nonlinear surface density 

To calculate a perturbed ring's surface density a, it is convenient to make the 'local' 
approximation, which assumes that the interparticle forces are exerted primarily by nearby 
ring particles that lie a small distance £ away from the perturbed particle, where i <^ a. 
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Because these perturbing particles reside on streamlines that are close by, one can ignore 
the streamlines' curvature that only occurs at great distances ^ £, and instead treat the 
streamlines as if they were straight wires having a linear mass density A (a) that is essentially 
constant^ about a streamline whose semimajor axis is a. Note, however, that perturbations 
of the ring also changes its surface density a, since these streamlines can be compressed or 
rarefied in the radial direction. To see this, let 6m = XSi be the total mass that resides 
along a segment of a streamline whose tangential length is 6i. If the ring were unperturbed, 
Sm = X6i = aoSaSi, where Sa is the streamline's radial width, so A = a^Sa is the streamline's 
linear density. If, however, the ring were perturbed, then 6m = XSi = a6r6i, i.e. X = ao6a = 
a5r where (j(a, 9) is the perturbed ring's surface density. Consequently, mass conservation 
requires cr(a,0) = aQ{a)/J where J{a,6) = 6r/6a dr/da in the limit that the streamline 
width 6a 0. Thus 

(To (a) 



where 



Or 



J 



=)1 



(39) 



(40) 



is a measure of the streamlines' degree of compression ( Borderies et al.lll985l ). Note that this 
quantity is sensitive to gradients in the streamlines' orbit elements, since 



da 



oa 



RL - im—R^ 1 e-'""^ 



d{ae) 
da 



+ ikae 



(41) 



where Rm = ae and the wavenumber k = —m{dCb/da) is the rate at which the streamline's 
longitude's of periapse varies with semimajor axis a. The above can also be written more 
compactly as 



^i(r]—muj) 



n'^ = qe- 

where q = \/(e + e')^ + {kae^ and tanr/ 
since d{ae)/da = e + de/dx = e + e', so 



kae 
e + e' 



J = 1 - qcos{(f) + r]) = l-^e [ge^('^+'')] = 3?e |l - ^[eaexp(i0)]| 



(42a) 
(42b) 

(43) 



where = m{6 — Og — uj). The angle r] indicates when the ring's response to perturbations 
takes the form of a spiral density wave. When the ring's responds is a tightly-wrapped wave, 
|A;ae| ^ |e + e'| and so \ri\ ~ 7r/2, while Ir^l ^ 1 when there ring's response is not wavelike. 



^Actually, A varies around the orbital ellipse by a small fractional amount that is of 0(e) <C 1, which 
may be neglected. 
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Note that q = \TZ'^\ is the nonhnearity parameter of iBorderies et aLl (jl985l ). who point 
out that adjacent streamhnes cross if their separation \6r\ exceeds \6a\ at any longitude 6, 
which would occur if g > 1. However previous studies have sho wn that a ring's internal 
forces tend to adjust the streamlines so as to avoid crossing (c.f., IShu et al.lll985l ). and we 
expect g < 1 here, too. 

The ring's response to perturbations is said to be linear when q <^ 1. In that case, the 
fractional variations in the ring's surface density, {a — (Jq)/(Jq ~ 7^J^e*™'*^^~^^\ are sinusoidal 
and of low amplitude. However, when the ring is nonlinear, q = \TZ'^\ is not small, and 
large variations in a (a, 9) can occur. Although those variations in a (a, 9) would be periodic 
in longitude 9, Eqn. fl39|) shows that they would not be sinusoidal. Lastly, note that Eqn. 
also implies that a nonlinear ring's surface density can drop by no more than 50% since 
< 1, which also occurs in the troughs of nonlinear spiral density waves (jShu et al. 



I7^' 



19851). 



If this study were instead interested in spiral density waves, we would then note that 
tightly wrapped spiral waves have wavenumbers \ka\ ^ 1 and amplitudes e(a) ~ constant. 
In this case, the second term in Eqn. fl^Tl) would dominate over the first, so 71'^ ~ iklZm- 
This is known as the tight- winding approximation, and it is valid when \kae\ » |e + e'| 
such that 7] ~ ±7r/2 in this case. If we were to insert this into Eqn. fl39|) and proceed further, 



we w ould then arrive at the theory for nonlinear spiral density waves (c.f. IBorderies et al. 
19861 ). which will in fact be considered in a foUowup study. 



However this effort is interested in the motion of particles orbiting near an inner Lindblad 
resonance, most of whom are on the non-wave side of the resonance. Also recall the results 
of Section 12.31 which suggests that the ring particles there will have similar longitudes of 
perihelia Co. Since uj varies slowly with radial distance x, the \kae\ term in the above is 
usually small in comparison to |e + e'| ~ |e'| term, since the streamline's eccentricities e(a) 



grow rapidly with proximity to the resonance. Thus TZ'^ ~ {de/dx)e 



qe 



vq 



e', ?7 ~ or 



TT, and a{a,9) = cxola) / J ^ ao/[l — e' cos m{9 — 9 s— uj)] in most (but not all) of the scenarios 
considered here. 

Finally, recall that the forced eccentricity of a single isolated ring particle is e{x) = ips/x 
(Eqn. I29ap . so the threshold for streamline crossing would be |7?.^| = \de/dx\ = ips/x^ = 
1. Solving for |x| = xnl then provides a rough estimate for the fractional distance from 
resonance where the ring's surface densit y variations will be no nlinear, and where the ring's 
internal forces should be significant (e.g.. IBorderies et al.lll982l ): 



Xnl 



3(m — e) 



(44) 
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This study is interested in Mimas' m = 2 ILR that hes near the B ring's outer edge, which 
has e = 1 and /™ = 1.50. The perturbing satelhte's mass is /i^ = 6.60 x 10^® Saturn masses, 
and its m = 2 resonance hes at = 117, 553.7km from Saturn's center, so tps = 3.30 x 10~®, 
xnl = 1-8 X 10~^, and Aajy^ = x^LO-r — 20km is the physical distance from resonance where 
nonhnear effects play a role in the ring's dynamics. We also note that the outer edge of 
the A ring is maintained by an m=7 ILR with the Janus/Epimetheus coorbital pair. That 
resonance has fj" = 9.4 9 , /x.^ = 3.32 x 10^^ (which is the mass of Janus, the larg er of the two 



coorbitals, iPorco et al.l (120071 )). and a mean semimajor axis = 136,773km (jPorco et al. 



19841 ) ■ so tps = 1-75 X 10~^, xnl = 4.2 x 10^^, and Aa^L — 6km. It is over these spatial 
scales where the ring's internal forces, such as self gravity, are expected to play an important 
role. 



2.4-3. ring gravity 

In the local approximation, a perturbing streamline is an infinitely long, straight wire 
having a linear mass density A(a') = (7o{a')6a', where cro(a') is the unperturbed ring's surface 
density, a' its semimajor axis, and 6a' is its radial width in semimajor axis space. The 
gravitational acceleration that streamline exerts on a nearby particle is Sag = 2GX{a')/d = 
2Gao{a')Sa' /d, where d = r' — r' is the separation between the perturbing streamline, whose 
planetocentric distance is r', and the perturbed ring particle at r. Since the streamlines' 
eccentricities e(a) are small, the streamlines are all nearly parallel, so the gravitational 
forces that they exert are radial. If the particle's semimajor axis is a, then that separation 
is d = r' - r = a' - a - 3?e{[/?„(a')e-™^'^' - /?„(a)]e^""(^-^=-^)} by Eqn. ([MD, where 
Au' = uj{a') — uj. This separation is then written more compactly as = a[x' — x)3?e[l — 
^gim(0-e,-<i)j ^ ^i^^i _ ^^i^Y - 5 COS 0), where = ra{Q -6^-00)+ (3 and 

^ ^ i?^(aOe"-^^' - Rm{a) ^ e{x')e-^^^^' - e{x) ^ 
a' — a x' — X 

where 

t /3 ^^{^') sin(mAco'') ^^^^ 
e(x') cos{mAuj') — e(x) 

and 6 = \A\ is real. The total acceleration that the entire ring exerts on the particle is then 

^ ^) _ /■ Sa - r°"' 2Gao{a')da' 

^ ' iring ^ Aij, a(x'-x)(l -5COS0) 

where a\n and Oout are the semimajor axes of those streamlines at the ring's innermost and 
outermost edges. 
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Next, Fourier-expand the longitude- dependent factor (1 — 5cos0) ^ in Eqn. P7|) . which 

112 

; = ^ + -[(1 - (5^)"^/^ - 1] cos(f)+ other terms, (48) 

1 — d)cos0 V 1 — 5^ d 

provided 6 < 1 such that streamhnes do not cross. Those other terms in the above are 
proportional to cos(ri0) = cos[nm(^^ — 6^ — u) + nf3], where the integer n takes all values of 
n > 2. Those other terms represent non- resonant forcings of the ring, and they are negligible 
since Rnm will be small. Inserting the surviving terms in Eqn. pHj) back into Eqn. P7|) then 
shows that the ring's gravity has the desired form of Eqn. fl7b|) : 

AIM ^ 2G r ^ , -°";'>;g_ (4ga) 

J a;,, (a;'-x)Vl-52 

where the r subscripts indicate that these are radial accelerations. 

For numerical work it will be convenient to put Eqns. (H9!) in a dimensionless form, 
which is obtained by dividing by \aV\ = 3(m — e)afi^. Also let s{x) = cr(x)/croo be the 
ring's fractional surface density, where the constant a^o is the ring's undisturbed surface 
density far away from the resonance. Also introduce the so-called normalized disk mass 
fid = iKTood^/Mp = TrGuoo/af^^, which is a dimensionless measure of the ring's unperturbed 
surface density that is quite convenient. With these definitions in hand, the dimensionless 
version of the above accelerations become 

A%. _ 2fid s{x')dx' 



|aX>| 37r(m - e) i^,.^ (x' - x)Vl - 5^ 

Am 

lis 
\aV\ 



and a^,{x) = — ^ = a^.e"*'"^ (50b) 



where = -f^^^ ^ "^'f , (50c) 

^"^^ ' 371 {m- e) A.^ e(x')e*™^'^' - e(x) ^ ' 

with the integration variable x' being the fractional distance from the nominal resonance 
that ranges over the ring's inner and outer boundaries Xin and Xout- 
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2. 4- 4- ihe hydrodynamic approximation 



Currently, there is some uncertainty in the rings hterature on how to handle the micro- 
scopic interactions that occur among ring particles. Those interactions are important, since 
they control the dynamical heating of a ring as well as the transport of angular momentum 
through the ring. For many years it had been thought that collisions were the most important 
particle-particle interact ion, and sophisticated models were developed to handle the effects of 
collisions. For instance, iBorderies et al.l (119831 ) solve the Boltzmann equation for the stress 
tensor that accounts for those interparti c le inte ractions, which yielded a formalism known as 
the particle-gas model. IBorderies et al.l (119851 ) also derive an alterna te stress tensor for an 
incompressible fluid-ring. We also note that the particle-jam model of iMosqueira fc Estrada 
( 120021 ) appears to be an intermediate case, since their ring is quasi-incompressible. 



However, we now know that self-gravi tating wakes also play an important role in stirring 
dense planetary rings (jPaisaka et al.ll200ll ). an d that the viscou s overstability can be another 
mechanism for viscous transport in the ring (jSalo et al. 2001). In light of the varied and 
uncertain ways in which ring particles interact on microscopic scales, we have elected to treat 
these interactions in the simplest way possible, via the hydrodynamic approximation, which 
assumes that the ring is a Newtonian fluid whose internal stresses manifest themselves via 
pressure an d viscosity. This hydr odynamic approach has been described as "qualitatively 
acceptable" (IBorderies et al. I ll982h . but we consider it to be the most reasonable and practical 
approach, given the current uncertainties. 



A Newtonian fluid is one where stress is proportional to strain. However, iLatter fc Ogilvie 



( 120061 ) use kinetic theory to show that the influence of w akes and other ins t abilit ies in a 'di- 
lute' planetary ring can in fact be non-Newtonian. But iLatter &: Ogilvid (120061 ) also show 
that the results of the kinetic theory do agree with hydrodynamics when the ring is dense and 
collisions are frequent, as is the outer B ring. It should also be noted that the many spiral 
density waves that are seen in Saturn's rings provides observational evidence that the hydro- 
dynamic treatment is in fact appropriate, since the damping of linear sp ir al density waves 



i s well described by simple hydrodynamic viscosity (ITiscareno et al.ll2007l ). ISremcevic et al. 



( 120081 ) also show that subtle variations seen in these waves' dispersion relation are probably 
due to hydrodynamic pressure. Both of these observations support the notion that the hy- 
drodynamic approximation is in fact appropriate when studying a planetary ring's large-scale 
motion, and justify its usage here. 
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2.4-5. particle collisions and pressure 

The following will account for the effects of particle collisions in a simple quantitative 
way, with coUisional effects being treated as if they were due to a pressure p(r, t) that tends to 
repel particles away from regions of higher density (where collisions will be more frequent and 
vigorous) and towards regions of lower density. To effect this, adopt a barotropic equation of 
state that assumes that the pressure in the ring can be expressed as the function p(r, t) = p{p) 
where p(r, t) is the ring's volume density that obeys = dp/ dp where c is the ring particles' 
dispersion velocity. Then the radial acceleration that an individual ring particle experiences 
due to coUisional pressure is 

1 dp Idpdp dp da 

^ p dr p dp dr p dr a dr 

since p = a /2h where h = c/VL is the ring's vertical half-thickness. Note, though, that this 
acceleration depends on a radial gradient, so this formula will be problematic at the edge of a 
sharp-edged ring. However that problem is avoided by first considering the linear momentum 
flux Gp that pressure effects also transmit through the ring. 

The linear momentum flux Gp is defined as the force-per-length (e.g., the two-dimensional 
pressure) that one streamline exerts on that exterior to it due to pressure effects. Thus 
5f = Gp6i is the radial force that one streamline exerts on a neighboring streamline segment 
whose tangential length is 6i. If those streamlines reside in the ring's interior (rather than 
at its edge), then the net force on a segment is A/ = 6f{r — 6r) — 6f{r) ~ ~{d6f /dr)6r = 
— {dGp/dr)5r5i where 6r is the segment's radial width. And since A/ = ap6m where 
6m = X6i = a6r6i, then the acceleration due to pressure is related to the linear momentum 
flux via 

. .1?^ . .L?^, ,52) 

a or (Jo oa 

noting that a = a^/J and that dr = J da. A comparison of Eqns. flSTl) and fl52|) shows that 
the ring's linear momentum flux due to pressure is 

Gp = c'a = (53) 

assuming that c is constant. 

Next, Fourier expand the longitude dependent factor in Eqn. (!53|) . which is Eqn. 
( |48|) with 6 ^ q and — > m{6 — 6^ — u) + r] where t] is from Eqn. fl42b|) . so 



Gp ~ , 



^l-g2 



1 + 



2(1 - a/1 - g^) j[m{9-9,-Cj)+r,] 



(54) 
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upon reverting to complex notation. Inserting this into Eqn. fl52l) shows that the acceleration 
due to pressure now has the desired form of Eqn. ([7]) where 



pr 



c2 d 



a^da y ^1 - g2 



A. 



pr 



2c^d_ 
(Jo da 



^ ^ ) i{r)—niCj) 



(55a) 
(55b) 



with the r subscript indicating that these are radial accelerations. The dimensionless ver- 



and 



where 



and 



are 






_ ^pr 


\aV\ 


oipj.{x) 


Am 

_ ^pr 


\aV\ 


ttpf,(x) 




3(r 


9p{.x) 


.(1- 



J2 



d 



3fm — t)s dx 



-m -imu) 



2 d 



12 



^ + ikagp 



q^ 



q' 



(56a) 

(56b) 
(56c) 
(56d) 



with c' = c/ a£l = h/a being the particles' dispersion velocity in units of their orbital velocity. 
The a™ quantity is introduced here for convenience, since it is a real quantity when the ring 
is nearly peri- or apo-aligned. Since \muj\ ~ or 180° in that case, qe^^ ~ e' in Eqn. (I56dp . 
while the terms proportional to kagp terms are negligible. Most (but not all) of the ring 
scenarios considered here are in this limit. 



Keep in mind that the accelerations in Eqns. are only valid for a ring particle 

that orbits in the ring's interior. Those equations do not apply to particles inhabiting the 
streamline at the ring's outer edge, since the pressure there is exerted only by the streamline 



that is orbiting just interior to it. In that case, 5f = Gp6i = 6map 



XttpSi is the force on 



a streamline segment having a tangential length 5i due to pressure, so the acceleration on 
the ring particles there is = Gp/X = Gp/cxQAa where Aa is the outermost streamline's 
semimajor axis width. Inserting Eqn. (15^ into then provides the coefficients for the 
acceleration due to pressure at the ring edge, whose dimensionless forms are 



^pri-^) 



and (x) 



AO 

\aV\ 

Am 
^pr 

\aV\ 



c'\a/Aa) 



3(m — e)^yi — 



-m -imuj 
^pr^ 



(57a) 
(57b) 
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a 



2 (i - ^/r^^ 

where cC{x) = ^ ^ ( ) e'"- (57c) 

This derivation also illustrates the utility of deriving the acceleration due to pressure from 
the linear momentum flux Gp, since this approach easily handles the discontinuous pressure 
drop that will occur are the ring's outer edge. 



2.4-6. ring viscosity 



The viscous acceleration a,y that a parcel of ring material experiences along one of the 
Cartesian axes will be denoted as a^i, where the i = 1,2,3 subscript indicates one of the 
X, y, z axes. That acceleration is 



- E- 



+ 



dxn dx. 



2 \- duk 



+ —CS^ ( — 



(5^ 



where p is the ring's volume d ensity, r? its shear viscosi t'^{§, ( is its bulk viscosity, and Ui is 
the velocity along the i^^ axis (ILandau fc Lifshita 119871 ). To transform this into cylindrical 
coordinates, replace xi — r so that the radial velocity Ui ^ = r and the tangential 
velocity M2 — >■ f e = t9 in the above, while the differentials dxi — Or and 8x2 rdO. Next, 
note that the ring's various quantities, such as its density or noncircular velocity, are all 
expected to vary rapidly in the radial j = 1 = k direction, but vary slowly in the tangential 
direction. Consequently, only the j = 1 = k terms need to be preserved in the above, so the 
components of the viscous acceleration simplify to 



d_ 

dr 



4 \ dVr 



pa^e 



d ( dvf, 



dr \ dr 



(59a) 
(59b) 



Note that Eqn. ( ]59b[) implies that a^e 7^ in a rigidly rotating disk having Vr = rO 
where 6 is constant. But this is unphysical, because it implies that a rigid rotator would 
also experience a viscous transport of angular momentum. However, this problem is easily 
fixed by replacing the dvg/dr in Eqn. (]59b|) with rdO/dr, which provides a more physical 



^Note that the shear viscosity 77 that is used only in Section [2.4.61 is distinct from the angle 77 of Eqn. 
(|42b| that appears elsewhere in this paper. 
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expre ssion that is also in agreeme nt with other treatments of viscous astrophysical disks 
(e.g., iLynden-Bell fc Pringldll974l ). Also replace the viscosities in the above with rj = Vgp 
and C = ^bPi where Vg is the kinematic shear viscosity and the kinematic bulk viscosity, and 
integrate the repaired version of Eqns. fISIJl) along the vertical direction so that j pdz — > a, 
which results in 



Ub a- 




dVr 

dr 



(60a) 
(60b) 



Note that the viscous acceleration also depends on radial gradients, which would be prob- 
lematic at a ring edge. But that difficulty is again avoided by considering the ring's radial 
flux of angular and linear momentum. 



angular momentum flux 

The viscous angular momentum flux Fy{r, 6) is the rate per-unit-length that one streamline 
transmits angular momentum to that orbiting just exterior to it via the ring's viscous friction. 
Thus 6t = F^6i is the torque that a segment of length 6i exerts on its exterior neighbor. If 
that segment orbits in the ring's interior, then it is also torqued by ring material orbiting just 
interior to it, so the net torque on that segment is At = 6t{r — 6r) — 6t{r) ~ —{d6t/dr)6r = 
— {dFu/dr)6i6r where 6r is the radial spacings between the adjacent streamlines. And since 
the net torque on this streamline segment is also At = ra^gSm where 6m = X6i = a6r6i is 
the segment's mass, this provides a relation between the viscous angular momentum flux F^, 
and the tangential acceleration a^e, 

1 dF, 1 OF, 

aue = ^ = (61) 

rcr or aao oa 



Comparing this to Eqn. fl60bp shows that the viscous angular momentum flux is 



_ z^sO'oa^ do 



^ _,^,r - ^ -- y^^. (62) 

The ring's angular velocity is 6 = Q + 6i Q[l + 2eexp(i0)] where = m{6 — Og — Co) 
[see Eqn. [38]), so 89 /da ~ dVt/da + 2(fi/a)(l - J) since 

d d 
1 - J = — ieae"^) - a— [eexp(i0)] (63) 
aa^ ^ aa 
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by Eqn. fHSj) when small terms of 0{e) are ignored. The derivative of Eqn. fllSal) provides 
dfl/da ^ —{3Q/2a)[l + (m — e)da'l/dx] where = A^/\aV\ = a° + is the axisymmetric 
part of the ring's gravity + pressure, so 36 /da ~ — (3fi/2a) [4J/3 — 1/3 + (m — e)da^/dx\. 
Inserting this into Eqn. fl62p then yields 

fl 1 r .dal' 
2i> sa^aVt { — — 1 — 3(m 



dx 



(64) 



to lowest order in the streamline's ecce ntricity e. This express ion is equivalent to the viscous 
angular momentum flux that is given in lBorderies et al.l (119821 ) when the gradient in the ring's 
radial acceleration, da'^/dx, is negligible. Next, Fourier-expand the longitude-dependent 
factors in F,y, which yields = [1 + 2g cos(0 + r])]/ (1 — g^)^/^+ other terms, with the 
expansion obtained from Eqn. fj48l) . so 



fu + 2q\U 



4/3 



1 -g2 



where fy is shorthand for 

U{x) = {l-q'r 



3/2 



1 + (m — e) 



dx 



-Q 
3^ 



(65) 



(66) 



and remembering to preserve only the real part of Eqn. (!65l) . 

The angular momentum luminosity through the ring is the integral about a streamline, 

Cu{x) = j Fydl = 37rUusaoa^n. (67) 



This is equivalent to that given in lBorderies et al.l (119821 ) when da^/dx = 0, which is probably 
true for most planetary rings. Note that C^, oc fu{x), and that fu<0 when the nonlinearity 
parameter q > q* where q* = \/3/2 ~ 0.866 (again, provided that da^/dx ~ 0; see Eqn. 
|66|1. This is the threshol d for the angular momentum flux reversal that was first described 



m 



Borderies et al.l (119821 ). When streamlines become so disturbed that q > q*, the viscous 



torque causes the streamline's angular momentum to flow inwards {Ci, < 0) rather than in 
the usual outwards direction. The ring particles orbiting in this disturbed region then spiral 
inwards due to this angular momentum loss, which also opens a gap in the ring. 



linear momentum flux 

The radial component of the ring's viscous acceleration a^r also transmits a flux of linear 
momentum where 

1 dG, 1 dG, 

aur = = (oo) 

a or (To oa 
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(see Section [2.4.51) . Comparison with Eqn. (I60ap shows that the ring's viscous hnear mo- 
mentum flux is 

'4 \ dvr f 4: \ ctq dvr 



G. = - ^-i^s + -^)-^ = - + ) Y^^- (69) 

The time derivative of Eqn. ([3SD provides the streamhne's radial velocity, which is Vj. = ri = 
-iuj^RmC'^ ~ -ieVLeae'^, so dvr/da ~ -ieVL{l - J) = -ief]ge^(<^+'') by Eqns. (I43|) and (l63|) . 
Inserting the real part into Eqn. (!69|) then yields Gu = — e + u^) flaoq sin(0 + r7) / J^. The 
Fourier expansion of the longitude-dependent factor in is sin(0 + ?7)/[l — g cos(0 + r^)]^ = 
2 sm{(j) + 7]) / [^y 1 — q'^{l + a/1 — q'^)] + other unimportant terms, so the viscous flux of linear 
momentum is 

~ , ' ^ (70) 

(i + vT^j 

upon switching back to the complex notation. 



acceleration due to viscosity 

Inserting the linear and angular momentum fluxes G^ and into Eqns. (pTl) and ( l68l) then 
provides the radial and tangential acceleration that are due to the ring's viscosity, which also 
have the same form as Eqns. ([Tj) with 



2a(jQ dx 
aao dx 



and 



AZ.(a) ^ 



2ien /4 



aao 



CToq ( fu - 

ox 



4/3 



aoqe 



i{r)—mu>) 



a - g2 1 + 



q^ 



(71a) 
(71b) 

(71c) 



where it is assumed that only (Tq, e, and uj might vary rapidly with distance x while the u^s 
are treated as constants. The dimensionless versions of the tangential accelerations are 



and 
where 

with 



a^e[x) 



a^0[x) 



9e{x) 



\aV\ 
A 



2(m — e)s dx 



\aV 



u9 —m —imu) 



sq f, 



m — e)s \ dx 
4/3 



dge , ., 
— h ikagt 



1 + v/r^ 



(72a) 
(72b) 
(72c) 

(72d) 
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where u'^ = Us/a'^Q is a dimensionless version of the ring's shear viscosity, and s{x) = 
a{x)/(Joo its fractional surface density. Similarly, the dimensionless version of the radial 



acceleration is 



Am 

and aZ{x) = ^ = ^^7^6"""" (73a) 
with fi'r(^) = — / C^3c) 



where z/^ = vi,/a?VL. Also note that the o^q and become real quantities when the ring is 
peri- or apo-aligned. 

Again, these accelerations are only valid for ring particles orbiting in the ring's interior, 
and do not apply to particles that inhabit the streamline at the ring's outer edge. In that case, 
5t = F,y6i = Smray = Xra,y6i is the viscous torque that is exerted on a streamline segment 
on the edge that has a tangential length Si, so the viscous acceleration on the ring particles 
there is a^, = F^/Xr ^ Fu/cTQaAa where Aa is the outermost streamline's semimajor axis 
width. Inserting Eqn. (1651) into a^, then provides the coefficients for the tangential viscous 
acceleration, 

«°.(^) = f% = :/^(7f) (74a) 
\aV\ 2{m — e) \Aa/ 

Am 

and aZix) = ^ = (74b) 

^ ::i^(/.-^^)(£,), (-) 

for particles orbiting in the ring's outermost streamline. Similarly, the radial component of 
the viscous acceleration is a^^ = G^/aoAa for particles in the outermost streamline, so the 
dimensionless version of this acceleration is 

aZ(^) = (75a) 

where ^Z(^) = ^iH^f^^"^ ( M . (,5^) 

3(m-e)v/r^(l + yr^j ^^""^ 



Finally, note that there are several known sources of visc osity in planetary rings: inter - 
particle collisions which have i/ oc cr in an optically thin ring (iGoldreich fc Tremaindll982l ). 
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and self-gravitating wakes that resuh in z/ oc cr^ fjPaisaka et al.ll200ll ). N body simulation s 



also indicate that the viscous overstability varies as v 



[see Fig. 13 of lSalo et al.ll200lf ). 



Alternatively, one could also have employed a more sophisticated viscous stress te nsor to 
account for the effects of collisions (iBorderies et al.lll983l : iMosqueira fc Estrada! l2002l ). How- 
ever, previous studies have shown that the s urface density near a perturbed ring-edge tends 
to be nearly constant (IBorderies et al.lll989l ). so the results obtained here are not expected 
to be particularly sensitive to any viscosity law. In light of this, a simple constant-viscosity 
law is employed here. Nonetheless, if a power-law z/ cx cr^ viscosity were instead preferred, 
then the formulas derived here may be adapted. 



24.1. drag 

Small ring particle s are also susceptible to dr ag forces, such as Poynting- Robertson 



(PR) and plasma drag (IGoldreich fc Tremaine 



19821). atmo spheric drag (IGoldreich fc Porco 



19871 ). and possibly the Yarkovsky effect (IRubincaml |2006| ) . Drag forces can cause particles 
to migrate radially, and are particularly effective at damping orbital eccentricities. Although 
large ring particles are relatively immune to drag forces, they can still experience the effects 
of drag indirectly by colliding with smaller drag-sensitive ring particles. 

Most drag forces vary with the particle's velocity relative to the local circular speed Av, 
so this work will assume that the acceleration on a ring particle due to drag has the generic 
form = — Cf^i^Av where Cd is a dimensionless drag coefficient whose value will depend on 
particulars of the unspecified drag force and the particle size. Since Av = ff + {r6 — aQ)6 = 



ief)eaf2e*'^ where = mOo + oOmt 



(from Eqns. [33] and [33]), this 



drag acceleration has the same form as Eqns. ([7]) whose radial and tangential components 

—leA^Q. The dimensionless versions of these drag 



are A 
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This drag force will damp the ring particle's eccentricity at a rate that may be obtained 
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by inserting the drag accelerations and A^q into the equation of motion fl2T]) . which yields 

rft ~ ~2* 

This corresponds to an e-damping timescale of 



:Cden. 



(77) 



de/dt 



P. 



orb 



(7J 



where Porb = 27^/^7 is the particle's orbit period. 



2.5. dimensionless equation of motion 

The equation for the particle's motion is Eqn. ([T]) with a = Rg + ap + a,^ + 3.^ being 
the acceleration due to the ring's gravity, pressure, viscosity, and a possible drag. Now that 
the ring's internal forces are suitably Fourier-decomposed (Eqns. |49l [55l and [7TH76|) . insert 
those as well as the anticipated solution for the particle's motion (Eqn. [33!) into its equation 
of motion (12T|) to obtain 

DTZm = -^T = -^7 + — - (79) 

where D{r) = — uj"^ is again the particle's frequency- distance from resonance. The 
dimensionless version of this equation is obtained by dividing by \aV\ and noting that TZm = 
ea exp{—imuj) and Um = e/t ~ ei7 near a resonance, so 

eede"™^ - 2iea^ + a'; + ^s = (80) 

where d{x) = D/T> is the dimensionless version of and the constant V's is the satellite's 
dimensionless forcing function, Eqn. ( I29bl) . The radial part of the dimensionless acceler- 
ation is the sum of the contributions due to ring gravity, pressure, viscosity, and drag, 
so aj" = A^/\aV\ = [a^ -|- -|- i{a^r + (^7r)\^~^^ i while the tangential acceleration 
is = v4™/|aP| = {o^q + a%)e-^'^ by Eqns. ^ ESHSZl ESHZHl and [76]). Then let 
= -|- a™, which is the sum of all the accelerations that are due to conservative forces 
(gravity plus pressure) while ipd = ~{p^T + "^X) + 2e(tt^ + aj^) becomes the sum of all the 
dissipative accelerations (viscosity plus drag), so that the equation of motion becomes 

{ede + cC- #d)e-*"" + = 0, (81) 

where ifjd is called the dissipative forcing function, in analogy with Eqn. (fT9|) . The factor 
d{x) = D/V in the above is obtained by inserting Eqns. (fTS!) into D to show that 

rf(x) = X - (3e- m)(a° +a° ), (82) 
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where x is understood to the fractional distance from the nominal resonance, which is where 
the resonance would be if the ring's internal accelerations where zero. Equation fl82l) thus 
shows how the axisymmetric part of the ring's radial acceleration will displace the resonance. 

The equation of motion (IHTl) is complex, so its real and imaginary parts provide two 
coupled equations for the streamline's forced orbital elements e(x) and uj{x). However those 
equations decouple when the dissipation is weak, \ipd\ ^ li^sl, which also results in a ring 
that is nearly peri- or apo-aligned such that \moj\ <^ 1. In that case, the accelerations a™ 
and ipd are all real and depend only on e(x), so the real and imaginary parts of Eqn. (ISTI) 
yield 



As one might expect, a numerical algorithm that uses the approximate decoupled equations 
of motion (183|) is will converge to a solution much faster than one that attempts to solve the 
exact coupled Eqn. (18T|) . But keep in mind that these approximate solutions to the equations 
of motion are only valid when \muj\ <C 1 and \kae\ <C \e'\. 

Lastly, we note that this problem also has a third unknown, the ring's surface density 
ao{a) or, equivalently, its fractional surface density s{x). To address this, a third equation 
for this quantity is derived below, which is obtained by requiring that all of the torques 
exerted on each streamline balance to zero. 



Two mechanisms transport angular momentum through the ring: the ring's viscosity, 
and the satellite's gravitational torque. The rate at which this transport occurs is the ring's 
angular momentum luminosity, and it has two parts, C{a) = C^ + Cs, where £y is the viscous 
angular momentum luminosity (Eqn. 1671). and Cg is the angular momentum luminosity that is 
due to the satellite's torque on the ring. Note that C must also be conserved, i.e., dC/da = 0, 
for if this were not the case, then streamlines would gain or lose angular momentum over 
time, which would also cause them to evolve radially since dC/da oc Ag. Consequently, 
static equilibrium requires C{a) to be a constant everywhere, which also provides that third 
equation for the ring's unknown surface density (To(ct)- 

To calculate the ring's angular momentum luminosity Cs that is due to the satellite 
torque, first consider the specific torque that the satellite exerts on a single ring particle, 
which is Ti(r, = r x (— V$s) ■ z = —d^g/dO = m(jf^{r) smm{9 — 9^). The time-averaged 



and uj 



eed + a™ + ~ 




(83a) 
(83b) 



2.6. angular momentum transport 
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torque on the particle is obtained by inserting r{t) = a + ri and 6{t) = Oq + + 6i into Ti, 
Taylor-expanding to first order in the small quantities ri and 6i (which are Eqns. [33] and 
inZj), and then time-averaging, which yields 

Ti ^ + [sinm(0-^,)+m^iCosm(^-0,)]^ = ^mi?^^^ sin(mt^), (84) 

where the brackets indicate a time-aver a ge oy er the forcing period 27T/\um\- An equivalent 



expression is also derived in iHahn et al.l (Il995l ). The time-averaged torque that the satellite 



exerts on an entire streamline is then STs = TiSm, where 6m = 2TTaoa6a is that streamline's 
mass and 6a its radial width. Another useful quantity is the satellite's radial torque density 
dTs/da = 6T/6a, so 

= mnaoaRm'^'^ sin(muj) . (85) 

oa 

Also note that the torque that a streamline exerts on the satellite is simply —6Ts, so the total 
gravitational torque that is exerted by all streamlines having semimajor axis (sma) interior 
to a is ^ 

Tr{a) = — f 6Ts = — mnaoaRm{a)'^^ sin {mCu {a )) da . (86) 

J sma<a J 

This torque is a second-order effect since it is the product of two small quantities — the 
satellite's forcing and the streamline's epicyclic response Rm- This explains why the 
ring-satellite torque was formally absent from the linearized equations of motion fll2ap . And 
since Eqn. (jHSj) is the torque that the ring interior to a exerts on the satellite, it is also the 
luminosity of angular momentum Cs that the ring material interior to a transmits gravita- 
tionally to the satellite, which becomes 

Cs{x) = = -em/™/i,/irf7(x)Mp(a^])2 (87a) 
where 7(x) = / s{x')e{x') sin{mLu{x'))dx' (87b) 



when \E'^ is replaced by Eqn. fl2^ . ctq = s{x)aoo, and R^ = e{x)a. 
The ring's total angular momentum luminosity is 

C{x) =C, + Cg = 3^sz/>rfMp(a^])2 - emf^^fisl^dlMpianf 



which is a constant everywhere when the ring is in static equilibrium. At sites far from the 
resonance, e and 7^0 while s and f,^ 1, so C = Su'^^dMp^aQY is the ring's angular 
momentum luminosity. This is also the rate at which the unmodeled part of the ring, which 
lies interior to the simulated region, delivers angular momentum to the simulated part of 
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the ring. Inserting this into Eqn. flHHl) and dividing by L then provides a third equation of 
motion for the third unknown, the ring's fractional surface density s[x): 

+ = 1 (89a) 

where iy[x) = sf^ (89b) 

and Ux) = -l!!^/^7(x) (89c) 

are the dimensionless angular momentum luminosities due to viscosity and the satellite's 
perturbations, in units of C Equation fl89al) can also be used to determine the location 
of the ring's edge, since the edge streamline at a; = Xedge is where the viscous torque is 
counterbalanced by the satellite's torque. Or, equivalently, the edge is where the ring's 
angular momentum luminosity £(a;edge) is entirely due to the satellite's gravitational torque. 
Consequently, the ring's edge is the streamline that satisfies isi^edge) = 1 and iui^edge) = 0. 

Lastly, we note that if the dissipation in the ring is weak, i.e., \muj\ <^ 1, then uj ^ 
ipd/Tnips (see Eqn. I83b[) . and the fractional angular momentum luminosity due to satellite 
perturbations becomes 

< J 

where tpd is the dissipative forcing function. Most of the ring scenarios considered below will 
be in this limit. 



3. Numerical simulations of a sharp-edged planetary ring 

The streamline model described above provides three equations of motion, the real and 
imaginary parts of Eqn. flHTj) and the torque-balance Eqn. fl89|) . which are to be solved 
for the system's three unknowns, the ring's eccentricity e(x), longitude of periapse a)(x), 
and the ring's fractional surface density s(x), where x is the radial coordinate in the ring. 
Note that this is a coupled set of nonlinear integro-differential equations, since the ring's 
gravitational acceleration requires integrating over the unknown e{x), Cj{x), and s{x) (see 
Eqns. I50!l . while the accelerations due to pressure and viscosity involve derivatives of these 
quantities (Eqns. [561 and [721— [73]) . The next subsection describes how these equations are 
solved numerically, with the streamline model then applied to the outer edge of Saturn's 
B ring. The following subsections then show how the simulated outcomes vary with the 
model's physical parameters, which are the ring's unperturbed surface density (Too, the ring 
particles' dispersion velocity c, the ring kinematic shear and bulk viscosities Vg and z/^, and 
the drag coefficient Cd- 



- 30 - 



3.1. numerical method 

To solve Eqns. flHTl) and fl89|) . treat the broad planetary ring as if it were composed of 
N discrete streamlines that are uniformly spaced with semimajor axes Oj = + [i — 1) Aa, 
where the ring index i ranges between 1 and N, and Aa is the rings' width in semimajor 
axis. The i^^ streamline's eccentricity is = e(aj), its longitude of periapse uji = uj{ai), and 
its fractional surface density Si = (yo{ai) / a^. This discretization then allows integrals like 
Eqns. (1501) to be replaced with sums over the Cj, cjj, and Sj, with finite differences used to 
calculate the derivatives that appear in the accelerations due to pressure and viscosity. This 
results in a coupled system of 3A^ nonlinear equations for the 3A^ unknown Cj, tDj, and Sj, 
which are then straightforward to solve for numerically. 

For example, the m = component of the ring's dimensionless gravitational acceleration, 



Eqn. flSOal) . is the sum = Yl^jti^ij where Gij = 2fidSjAa/3iT{m — e){xj — 1 — 6fj 

is the gravitational acceleration that streamline j exerts on a particle in streamline i, where 
fid is the ring's normalized mass (see Section [2.4.31] , and with 6ij calculated via a similarly 
quantized version of Eqn. (HSj) . Note that Ga = here, because in the local approximation, 
a straight wire-thin streamline exerts no net gravitational force on the particles inhabiting 
that streamline. A similar strategy is also used to calculate the dimensionless m > 1 part of 
the ring's gravity, a^(xj); see Eqn. (]50bp . Also note that the accelerations due to pressure 
and viscosity, Eqns. (13^ and (1721— [751). require derivatives of e, lj, and s, which are c alculated 



numerically using a three-point Lagrangian interpolation scheme ( lHildebrandlll956l ). Lastly, 
keep in mind that Eqns. fl57|) and flTiHTSl) are the accelerations on the N^'^ streamline at the 
ring's edge that is exerted by the adjacent — 1 streamline, so the right hand sides of those 
equations are evaluated at xn^i = (otv-i — ar)/ar. 

The number of streamlines A^, as well as their radial widths Aa, are both chosen so that 
the streamline model can readily resolve the disturbances that Mimas excites at the outer 
edges of the B ring. The following subsections will show that these disturbances are usually 
confined to the rings' outermost 50 or so km, so most calculations use A^ = 300 streamlines 
that have radial widths of Aa = 0.5 km so that the total radial extent of the simulated region 
w = NAa is usually about w = 150 km. For most ring models, this width is broad enough 
to show that the simulated ring's innermost part furthest from the resonance has in fact 
adopted the single particle solution, Eqn. fl29ap . This indicates that the ring's internal forces 
there are negligible, and that a wider a ring need not be simulated. However, the following 
simulations will also show that the self-gravity in a more massive ring is able to transmit 
the satellite's disturbance to greater distances inwards of the resonance. To account for the 
greater reach of these ring's internal forces, the streamlines' radial widths are increased to 
Aa = 1 or 2 km so that the total radial extent of the simulated region becomes w = 300 to 
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600 km in these more massive planetary rings. 

The streamhne model developed here, called NLsgvp, is written in IDL, and it solves the 
coupled set of 3A^ nonlinear equations for the streamlines' 3A^ orbit elements and surface den- 
sities. Although IDL does supply routines for solving coupled systems of nonlinear equations, 
they are not particularly robust, and sometimes fail to find a satisfactory solution to these 
equations. So the MPFIT algorithm is instead used here, which is an IDL procedure written 
by Craig Markwardt that is available at ||http:// cow. physics. wise. edu/~craigm/idl/fitting.html[ 
MPFIT is a parameter search algorithm, and it was formally designed to search parameter 
space for the purpose of fitting a parametrized model to a dataset. However, the act of seek- 
ing solutions to a coupled set of nonlinear equations is conceptually very similar to modeling 
data, since we wish to find the set of 3A^ 'parameters' Cj, u)j, and Sj, that best satisfy the 
3N equations fIST]) and fl5^ . whose right-hand sides are 'data'. MPFIT is quite efficient and 
very well-suited for solving this type of problem. 

MPFIT must be initialized by first providing it with a trial solution. That initial guess 
is formed from the single particle solution, Eqn. fl29ap . which is also the solution to the 
equation of motion fl83al) at sites far from the resonance where the conservative part of 
ring's internal acceleration, a™, is negligible. Note, however, that Eqn. fl29ap diverges at the 
resonance where d{x) ~ x — * 0, whereas e{x) is expected to stay finite inside the nonlinear 
region where < xnl (see Eqn. HI]) due to the ring's internal forces. So to qualitatively 
mimic this effect, replace x in Eqn. fl29al) with x —e^x'^^ x%^ such that 

etriniix) = — , (91) 

where ips is the satellite's dimensionless forcing function, Eqn. (]29b|) . This is the trial solution 
that is used to initialize MPFIT, which also adopts s{x) = 1 and uj{x) = 0. 

All of the calculations described below use double-precision arithmetic. To assess the 
accuracy of the model results, we insert all numerical solutions obtained by NLsgvp back 
into the equation of motion flSTl) and divide by ips, with those solutions also being inserted 
into the torque-balance equation fl89al) subtracted by 1. The residuals on the right-hand side 
of those equations are the solutions' fractional errors, which for the models described below 
are all smaller than 4 x 10^^^. 
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3.2. simulations of the outer B ring 

3.2.1. variations with surface density 



Figure d] illustrates how the ring's epicyclic amplitude Rm{x) = ae{x) varies with a^o, 
which is the ring's undisturbed surface density far from the edge. These curves are solutions 
to Eqns. flHTj) and fl89|) . and were obtained using the numerical method described in Section 
13.11 The main point of this figure is to show that the epicyclic amplitude Rm at the ring 
edge is larger when the ring has a lower mass. Note that the ring's response approaches 
the single particle result, Eqn. (!28l) . as cXoo — > 0. All of these simulations adopt viscosities 
of z/.s = = 50 cm^/s e c, which is compa rable to the viscosity measured in Saturn's A ring 
( iTiscareno et al.l 120071 : iPorco et al.l 120071 ). There is no drag in these simulations (C^ = 0), 
and the ring's outer edge is chosen to be at the nominal resonance at a: = = a — a^. The ring 
particles' dispersion velocity c is also chosen such t hat the ring is marginally gravitationally 
stable. A gravitationally stable ring has a. Q > 1 (jToomrelll964l ). where Q ^ cQ/nGcXoo = 
d I fid- However, Saturn's main A and B rings also exhibit brightness asymmetries that are 
believed to b e due to t he presence of self-gravitating wakes, but those wakes only occur 
when Q < 2 (js^ll992h . Consequently, a higher-mass but marginally stable ring also has 
a larger dispersion velocity, and the simulations reported in Fig. [T] have c chosen so that 
Q = 2. The dotted horizontal line also indicates the B ring's observed epicyclic amplitude, 
i?obs ^ 45 km, which was rn easured in Cassini spacecraft observations that were acquired in 
2005 ( ISpitale &: Pored l2006l ). Figure [T] tentatively suggests that the outer edge of the B ring 
might have a surface density as low as cToo ~ 10 gm/cm^. 

Figure [2] plots the nonlinearity parameter q versus distance from resonance for the 
simulations described in Fig. [H In all of these simulations, \kae\ -C |e'|, so the nonlinearity 
parameter q ~ |e'| is also the ring's eccentricity gradient. Figure [2] shows that the disturbed 
edges of lighter rings tend to be more nonlinear than heavier rings. The ring's epicyclic 
amplitude Rm is also proportional to the satellite's mass, so increasing the perturber's mass 
also makes the ring more nonlinear. 



3.2.2. variations in the ring-edge location 



Section 12.31 notes that the outer edge of the B ring could lie as far as 24km exterior to 
Mimas' nominal m = 2 ILR. This displacement was first measured by iPorco et al.l (119841 ) 
from Voyager observations, but with a large uncertainty. Ho wever, a prelimi n ary a nalysis 
of Cassini observations of the B ring, which is described in ISpitale fc Pored (120061 ). does 
indicate that the ring-edge lies about 20 or so km beyond the resonance. 
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Figure [3] illustrates the consequences of displacing the B ring's edge radially outwards; 
the three curves show the ring's epicyclic amplitude Rm{a) for simulated rings whose outer- 
most semimajor axis lies a distance Aa = 0, 10, and 25 km exterior to the nominal resonance 
location. Figure [3] also shows that the ring's epicyclic amplitude Rm{,o) grows linearly as 
the edge is approached, even in the portion of the ring edge that lies exterior the resonance. 
Consequently, letting the ring edge lie farther beyond the resonance will result in a larger 
epicyclic amplitude. However, Section 13.2.11 shows that the increase in the epicyclic ampli- 
tude can be offset by increasing the ring's surface density a^o- So Fig. [3] assigns a larger 
surface density to rings whose outer edges are displaced farther outwards, but in a manner 
that keeps the epicyclic amplitude at the edge comparable to the observed value, Rohs — 45 
km, which is dotted line in Fig. [31 The rightmost curve in Fig. [3], which agrees with the B 
ring's observed epicyclic amplitude as well as the likely upper limit on the edge's location, 
now shows that the B ring could have a much higher surface density of a^o ~ 280 gm/cm^. 



3.2.3. variations with dispersion velocity 

FigureHK shows that the ring's epicyclic amplitude Rm is quite insensitive to the particles 
dispersion velocity c. Those simulations are for a B ring that has a surface density fixed 
at (Too = 30 gm/cm^ and a stability parameter that ranges over < Q < 12, which is 
a convenient proxy for c since Q ^ cQ/irGaoo. These Q values correspond to dispersion 
velocities of c < 4.9 mm/sec and vertical scale heights h = c/Q of /i < 32 m. Note that 
the black Q = curve, which represents a pressureless ring, is indistinct from the other 
Q < 5 models. Figures Hb and plot the simulated rings' fractional surface densities 
s(a) and longitudes of periapse oj{a) versus semimajor axis a, which is able to distinguish 
the dynamically hotter Q ^ 10 models from the cooler Q < 5 models. Note the large 
variations in s{a) and uj that is seen at the outer edges seen of the hotter rings that have 
Q > 10; those variations are due to the larger pressure drop that occurs across the outermost 
streamline (e.g., Eqns. EZD- Those surface density excesses seen at the rings' edges are also 



reminiscent of those occurr ing in models of narrow ringlets (e.g., IChiang fc GoldreichI 12000 



Mosqueira fc Estrada! |2002| ) . but the variations seen in Figs. H] are poorly resolved due to 
the models' radial sampling of Aa = 0.5 km. Nonetheless, Figs. H] do show that pressure 
effects, if important at all, are confined to the ring's outermost Aa ^0.5 km. Since this 
study is interested in the perturbed ring's state over a much broader radial scale, Aa ~ 40 
km according to B ring models of Fig. HI the following will for simplicity keep Q = 2 for 
most of the ring scenarios that are considered below. 

Lastly, note that oj{a) > in the hotter Q = 12 model (Fig.HJ:). Equation (l85l) indicates 
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that the torque that the sateUite exerts on this model ring is positive. But that torque has 
the wrong sign, since a sateUite must exert a negative torque on the ring if it is to be confined 
interior to an ILR. This inabihty of the satelhte to confine a viscous and dynamically hot 
ring is another reason why hot ring models need not be considered further. The confinement 
of the B ring by Mimas' gravitational torque is also considered in greater detail below. 

3.2.4- '^^'^fl' viscosity, and the torque-balance problem 

Figure [5] shows how the simulated outcomes depend on the ring's viscosity. These models 
are similar to the B ring simulations shown in Fig. HJ except that the rings' viscosities take 
values of I's = = 5,50, and 500 cm^/sec while Q = 2. Note that the rings' epicyclic 
amplitudes Rm = ea are all insensitive to the viscosity = (Fig. [5^), as are the rings' 
fractional surface densities s (Fig. [5]d), and their angular momentum luminosities and ig 
(Fig. 01). This is due to the ring being nearly peri- or apo-aligned, i.e., \uj\ ~ \ipd/fn'^s\ ^ 1 
(see Eqn. I83bl) in all of these simulations. Since the dissipative forcing function ipa is small 
compared to the other terms in the equation of motion, it has little effect on e (see Eqns IHT] 
and [83]). However, Fig. [St demonstrates that the streamlines' longitudes of periapse u oc ipd 
are proportional to the ring viscosity, since 

i^d = -aZ + 2e«re = ^u^ i^s + fi^b) (92) 
where the factor / is of order 1 (see Eqns. [72]473|) . 

Figure [5li also shows that the ring's angular momentum luminosity is-, which is due 
to the satellite's gravitational torque on the ring, is small, i.e., \C.s\ 1- This is very 
problematic, because Section showed that the ring's outer edge should be the site where 
the viscous torque is counterbalanced by the satellite's torque, which is also where — > and 
£s ^ 1. However, this torque-balance requirement is not satisfied by any of the simulations 
described in Figs. [THSl Although the simulations described above provide useful illustrations 
of how a ring's simulated outcome depend on its physical properties cXoo, z^, and c, they are 
all unphysical since they do not achieve a torque-balance at the ring's outer edge. However, 
Section 13.2.51 does explore an alternate ring scenario that does in fact achieve the desired 
torque balance. 

But first, a final comment on Fig. Eb, which shows that the simulated ring's surface 
density has a ~ 20% excess in its outermost ~ 20km. This is a due the conservation of 
the ring's angular momentum luminosity £ ~ oc sfu (see Eqn. |88] and note that Cg 
is negligible in these simulations) where the function varies approximately as 1 — 4g^/3 
(see Eqn. [661) where the nonlinearity parameter is q ~ \de/dx\, which is also plotted in Fig. 
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|5ll. Because the ring' eccentricity gradient gets large near the resonance, the quantity 
is diminished there, but that must be offset by an increase in the ring's fractional surface 
density s = cq (a) /cr oo in order to conserve C oc sf^. 



3.2.5. possible B ring solution 



Section [3.2.41 notes that all of the ring simulations considered thus far fail to balance the 
ring's viscous torque against the satellite's gravitational torque, so the proximity of the B 
ring's outer edge near Mimas' m = 2 ILR is not yet accounted for. Equation ( IHTI) shows that 
the satellite's torque on the ring is proportional to 7(x) = se sin(ma;)c/x', which is the 
product of the streamline's forced eccentricity e(x) and the angle by which that streamline 
lags behind the satelhte's longitude, uj{x). That lag angle depends on the ring's shear and 
bulk viscosity via smmuj oc ipd oc + fut where / ~ C(l) (Section l3.2.4p . so the satellite's 
total torque on the ring is controlled by the sum of these viscosities. Although an increase in 
the ring's shear viscosity does increase the satellite's torque, this is of no help here because 
the viscous torque (Eqn. [67j) would increase by the same factor, and the torque mismatch 
would still persist. However, Eqn. flHTl) indicates that an increase in the ring's bulk viscosity 
z/fe can increase the satellite's torque without altering the viscous torque. 

But first, a comment on measurements of the viscosity of planetary rings. The shear 
viscosity Ug is, as the name suggests, a measure of the friction that results from the ring's 
shearing motions. On the other hand, the ring's bulk viscosity is the friction that occurs 
due to the ring's compressive or decompressive motions, which explains why Uf, enters into the 
equations of motion when there is a gradient in the ring's radial velocity (Eqn. [S9a|l . One way 
to infer and u^, is to examine spacecraft observations of spiral densit y waves; these waves are 



damped by visc ous effects over the radial scale Ar^ oc (z/^ + |z/ft) (iGoldreich fc Tremaine 



1978 



Shu 



19841). so a me a surem ent of Ar^ provides an estimate of the combined viscosities 



Tiscareno et al.l (120071 ) use Cassini observations of spiral de nsity waves in S a turn's 
A ring to infer viscosities of 30 < z/ < 300 cm^/sec there. However, iTiscareno et al.l (120071 ) 



dropped the bulk viscosity Ub from their Eqn. (7), so the viscosities quoted there should 
instead be interpreted as the combination u z/^ + |z/f,. Similarly, E gns. ([73ll an d (1751) 
indicates that our simulations will be sensitive to the combination I'g+^i'b- iPorco et al. (120071 ) 
also use the widths of the Keeler and Encke gaps in Saturn's A ring to infer viscosities of 20 < 
z/ < 87 cm^/sec, but again these viscosities are likely measurements of some combination of 
Ug and z/fe. The upshot is that current ring observations only constrain the sum of Pg and Ub, 
and that Ug and Ub are not known individually. 



With this in mind. Fig. M shows the results of three B ring simulations whose viscosities 
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sum to Us + i^b = 100 cm^/sec (which is roughly the A ring's total viscosity), while their 
ratios obey Ub/^s = 0, 1, 1000, and 7873. That last ratio was chosen so that the viscous and 
satellite torques do indeed balance at the ring's outer edge, which is where s ^ while 
i,^ —>■ and £s — > 1 (black curves in Fig. and d). These simulations' other parameters 
are similar to that adopted previously in Fig. [5l with cxoo = 30 gm/cm^, Q = 2, Cd = 0, and 
with the ring's outer edge stationed at the nominal resonance position at a = a^. 

But note that the B ring's mean outer edge actually lies about 24km exterior to the 
nominal resonance, and that an outwards shift in the resonance position would also result in 
a larger epicyclic amplitude Rm (Section 13.2.21) . But that can compensated for with a larger 
surface density (Too, which reduces Rm (Section l3.2.ip . Figure [7] illustrates one possible B 
ring solution that does achieve a balance of the viscous and satellite torques at the ring's 
outer edge at a — a,. = 24 km, which is where s while iu and 1, and also 

has the observed epicyclic amplitude of Rm — 45 km at the ring's edge. This simulation has 
(Too = 226 gm/cm^, Q = 2, Cd = 0, a. total kinematic viscosity of Ug + = 100 cm^/sec, a 
shear viscosity of Us = 0.00603 cm^/sec, and a viscosity ratio of Ub/us — 1.6 x 10^. Note that 
these simulations depend only on the viscosity ratio Vb/vs. For instance, when the simulation 
of Fig. [7| is executed again using shear and bulk viscosities that are xlO larger, the same 
Figs. b, and d are obtained, while the longitude of periapse uj (Fig. [Tt) is larger by x 10, 
as expected (e.g., Eqn. I83bl) . 

Figure Wp shows that the simulated ring is also very nonlinear, with q = 0.868 at its 
outer edge. This is actually just slightly larger the expected maximum possible value of 
q* ~ 0.866, but keep in mind that this q* threshold is approximate; see just below Eqn. (!67|) . 
This large nonlinear parameter also results in significant longitudinal variations in the ring's 
surface density a. These variations are demonstrated in Fig. [8], which plots radial profiles 
of the ring's surface density <j{r,9) along the satellite's longitude 9 = 6^, as well as along 
longitude 6 = Og ± 45° and 6 = 6s i: 90°. Those curves are calculated via 

a{a,e) s{a) ^^^^^ 



(Too 1 - g COS(0 + 7]) 

Ar{a,6) = r — ~ 0^(3; — e cos </>) (93b) 
and (j){a, 6) = m{6 — 6s — Cj) = (p* + ^^^(acdge) — ^{cl)] (93c) 

which may be obtained from Eqns. (!34|) . (139|) . and fH3|) . where Ar(a, ^) = r — is the radial 
distance of streamline a at longitude 6 from the nominal resonance a^, and a}(aedge) is the 
longitude of periapse at the ring's outer edge. Setting the angle = m[6~6s — cu(aedge)] = 
generates a radial profile along the ring-edge's longitude of periapse, while setting 0* = ±180° 
results in a surface-density profile towards the ring-edge's apoapse. Figure [8] shows that the 
periapse profile has a large surface density excess at the ring's outer edge, which is due to the 
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satellite's perturbation having shoved ring material inwards and compressing the streamhnes 
there. Conversely, the outer edge of the apoapse profile shows a broad low-density shoulder, 
which is due to the streamlines being more distended there. Since the |a;| and |?7| are all small 
in this model, the intermediate (j)* = ±90° profile, which is along longitude 6 = Og ± 45°, is 
also where J ~ 1 and r ~ a, so the surface density along this longitude equals its so-called 
undisturbed surface density o"o(a) that the ring has adopted in order in order to achieve a 
torque balance. 

So to summarize the results of this model of a viscous gravitating B ring: if the B ring is 
indeed a viscous Newtonian fluid, then its edge near Mimas' m = 2 ILR can only be accounted 
for when the ratio of the ring's bulk/shear viscosities at the ring's edge takes extreme values 
of lyb/^s ^ 10^. The ring's bulk viscosity Uh is a measure of the friction that results when 
the ring is compressed radially by the satellite's perturbation. When that is large enough, 
the shepherding torque exerted by the satellite is then strong enough to counterbalance the 
ring's viscous torque, thereby confining the ring in the vicinity of the resonance. However, it 
was a surprise to find that a more conventional model, one having a z/f, that was comparable 
or less than Ug, failed by a wide margin to balance the viscous and satellite torques. Also 
keep in mind that the ring viscosity u t hat is inferred from studies of planetary rings (e.g.. 



Tiscareno et al.l 120071 : iPorco et al.l 120071 ) is actually a linear combination of Ug and z/f,. So 
if the preceding scenario is correct, then the ring-edge's viscosity is dominated by its bulk 
viscosity Uf,, while the ring's shear viscosity Ug is negligible in comparison. If this finding is 
correct, then the B ring can shear freely with little dissipation, while its compressed radial 
motions are very dissipative. 

This is an unexpected finding. Nonetheless, the viscosity requirement Ub/ug > 10^ might 
be satisfied if the B ring's compressed regions are so densely packed that the ring particles 
there are 'shoulder to shoulder', with little voidspace between them. In this case, th e ring 's 



volumetric density would be nearly incompressible, as envisioned by lBorderies et al.l (119851 ). 
This volumetric incompressibility means that as the satellite attempts to compress the ring's 
streamlines further in the radial direction, ring particles must rise vertically and roll or tumble 
over each other as they are shoved inwards. This could be a very lossy process, one that 
might satisfy I'b/^s ^ 10"^, since the particles at the ring-edge would effectively experience 
2m 'avalanches' during each orbit. 

However it is unknown whether the compressed regions in a confined ring-edge are in 
fact close-packed. And the suggestion that tumbling close-packed ring particles might satisfy 
^b/^s ^ 10*^ is at this stage mere speculation. Indeed, the reviewer of this paper suggests that 
the non- Newtonian properties of self-gravitating wakes might inste ad play an import a nt rol e 



here. However, the relevant physics, such as the kinetic theory of iLatter &: Qgilvid (120061 ). 
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or the results of Nbody simulations (jPaisaka et al.l l200ll : ISalo et al.l l200ll ). are not easily 
adapted to a semi-anlytic treatment like ours, and so the consequences of non-Newtonian 
behavior are not explored here. In light of these uncertainties, we also consider an alternate 
ring confinement mechanism below, to demonstrate that there may be more than one way 
for a satellite to resonantly confine a sharp-edged planetary ring. 



3.2.6. drag in a planetary ring 

Section 12.4.71 notes that small ring particles are also susceptible to drag forces, such as 
PR, plasma, and atmospheric drag, and possibly the Yarkovsky effect. Although a drag force 
tends to have the greatest influence among smaller ring particles, the smaller particles can 
still communicate the effects of this drag to the larger particles via collisions. This is due to 
the enhanced eccentricity damping that a smaller particle experiences (Eqn. [77|) : its lower e 
will then favor collisions with the larger, more eccentric particles. The cumulative effect of 
these collisions would then resemble a drag force that also acts on the larger particles, too. 
But this scenario only works if the smaller ring particles are also sufficiently abundant, i.e., 
if the ring particles' size distribution is sufficiently steep. Obviously, the outcome of this 
scenario will depend on the particular drag force that is operative in the ring, as well as the 
particles' size distribution. And this initial study of drag is not prepared to deal with these 
uncertain details in a rigorous, quantitative way. Nonetheless, the effects of drag can still be 
assessed qualitatively by using the simple generic drag acceleration that is described Section 
12.4. 7[ and by assuming that there is a single effective drag parameter Cd that adequately 
describes how this drag effects the entire ensemble of ring particles. 

When the drag is the dominant source of dissipation, a particle's dissipation function 
becomes 

i^d = ^drag = -a7r + '^^a^ = (94) 

(see Sections 12.4.71 — 12. 5p . This drag force also causes the streamlines' longitudes of perihelia 
to trail behind the satellite's longitude 9s by an angle a) ~ ipdra.g/^4's = —^Cde/mf^^s (from 
Eqns. I29bl and I83bp . The torque Tg that the satellite exerts on the ring is — Ix Eqn. fl87ap . 
so Tg oc J^j^ esesmmojdx' and sgn(Ts) = — e. Consequently, this drag force also enables a 
shepherding torque that tends to drive ring particles radially away from the satellite's orbit, 
since < at an e = +1 ILR and > at an e = — 1 OLR. And if the shepherding torque 
is strong enough, it can counterbalance the viscous torque and maintain the ring's edge near 
the resonancqfl. This torque-balance is also known as resonance trapping, and it can occur in 



■^Of course, drag forces such as plasma drag, Yarkovsky, etc., can also exert an axisymmetric to = 
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a wide variety of disk-perturber syst ems that also have s ome dissipation, such as dust trapped 
at a planet's OLR due to PR drag (IRoques et al.lll994l ). and the ca pture of planetesimals at 
a protoplanet's OLR due to the solar n ebula's aerodynamic drag (IWeidenschilling fc Davis 
19851 : lMalhotralE993l : Iflahn et al.l[l995h . 



Ring confinement due to drag is illustrated in Fig. [9l which is for a B ring that has 
a surface density a^o = 265 gm/cm^ far from the resonance, a kinematic shear viscosity of 
Vg = 46 cm^/sec, and a drag parameter = 1-0 x 10~^. These parameters were chosen 
so that the epicyclic amplitude at the ring's outer edge has the observed value of Rm = 45 
km, and that the ring's viscous torque balances the satellite's torque at a distance of 24 
km beyond the nominal resonance. Keep in mind that Fig. [9] represents just one possible 
solution for the B ring, since other judiciously chosen combinations of ring parameters will 
also lead to a torque balance at 24 km beyond the nominal resonance. 

Figure ^ also shows that when drag is the dominant source of dissipation, then a large 
longitudinal offset becomes possible, with u ~ —23° at the ring's outer edge in this example. 
But keep in mirid that u oc —Cd, so stronger drag would result in a larger u. Interestingly, 
Spit ale fc Pored (120061 ) measured an offset of ~ —28° in Cassini observations of the B 



ring's orientation, yet Voyager observations revealed no significant offset (jPorco et al.lll984j ). 



Figure [TO] uses Eqns. fl93|) to calculate the ring's relative surface density a{r,9)/aoo as 
a function of radial distance r along selected longitudes 9 = Og (which is along the ring's 
longitude of periapse), 9s ± 90° (towards the ring's apoapse), and along the intermediate 
longitude 9s ± 45°. Again, the model predicts a large surface density enhancement at the 
ring-edge's periapse, as well as a low surface density shoulder at the ring-edge's apoapse 

Evidently, one of this model's main predictions is that the ring-edge's surface density 
should increase by ~ 50% at periapse; see Figs. M and [101 Presumably, this periapse en- 
hancement would have been seen by the Cassini spa cecraft, if not by Vo yager. One possible 
explanation for this non-detection can be found in iPorco et al.l (120081 ) . which describes a 
sophisticated photometric model of a swarm of ring particles that are illuminated by the 
Sun and Saturn, and also imaged by a nearby spacecraft. That model predicts that the 
optical surface brightness I/F of Saturn's main rings should saturate when the ring's op- 
tical dept h exceeds about 0.3 . Since the optical depth of the outer B ring exceeds that 
threshold (IColwell et al.l 120071 ). this I/F saturation could account for the absence of any 
detection of the expected periapse enhancement of the ring-edge's surface density. However, 
the ring should be less optically thick when observed at longer wavelengths, so the antic- 



torque on a ring particle that would cause it to migrate radially. In that case, a sharp ring-edge would 
instead indicate a balance between the satellite's shepherding torque and the m = part of the drag torque. 
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ipated periapse enhancement might be detectable when observing the ring-edge during a 
radio occuhation. 

The main purpose of this Section is to illustrate how other forms of ring dissipation 
might enable a satellite to confine a ring at a Lindblad resonance. For instance, drag 
forces might play a role in the confinement of the inner Uranian rings, which do orbit close 
enough to that planet to experience a drag with that planet's extended upper atmosphere 
( iGoldreich fc Pored 119871 ) . However, it should be noted that the differential particle size 
distributi on in Saturn's B ring vari es as dN{s) oc s~'', where s is the particle radius and 
g ~ 2.75 (IFrench fc Nicholson! l2000l ). Consequently, the B ring's differential mass distribu- 
tion varies as dM{s) oc s3-q~o.25^ which is weighted towards the larger particles, so the B 
ring is probably immune to the effects of drag forces. 



3.2.7. the satellite's torque 

As was noted above in Section 13.2.41 all of the B ring models that are described in 
Sections [32IIH3i23] and Figs. [THS] failed to find a solution that balances the satellite's torque 
on the ring against the ring's own viscous torque. This was rather surprising since most 
of those models were constructed (or so we thought) so that the satellite's torque would 
counter-balance the viscous torque. The maximum torque that a satellite can exert on 
material orbiting at its m**^ Lindblad resonance is 



\V\ ~ 3(1 -e/m)^'^^"^^ 



Tgt = - — y-"-^ = ^;\ fi,M,{any (95) 



( iGoldreich fc Tremaind Il978l . Il982l ). The viscous torque is the ring's viscous angular mo- 
mentum luminosity, C = Si^Ug/ a^fl) fj,dMp{afl)^ (see just below Eqn. [88]) . so one might expect 
a torque-balance to be possible when |Tgt| > Tiy, or when 



<fri^s? 

9(1 - e/m) 



a^n, (96) 



which evaluates to < 145 cm^/sec for Mimas' m = 2 ILR. Thus a torque-balance at the B 
ring's outer edge should seemingly be possible if the B ring's shear viscosity were comparable 
or less than that inferred for the A ring. However, the following review of the derivation of 
the Goldreich-Tremaine torque formula, Eqn. ( !95l) . will show why it can overestimate by a 
large margin the torque that Mimas exerts at a ring's outer edge. 

The easiest way to derive Eqn. (!95l) is to consider a drag-dor ninated ring whose dissi 



pation function is Eqn. (1941) : that scenario was also considered in iMeyer-Vernet fc Sicardy 



-41 - 



( ll987l ). The real and imaginary parts of a ring particle's equation of motion flHTj) are 



smmu; = — and cosmu; = , (97) 

where it is assumed that is real. Squaring and summing Eqns. fl97j) also shows that 
{eed + a™)^ + i^d — V's^ where ^/'d is Eqn. fl9^ and the is the acceleration on a particle 
that is due to gravity plus pressure. If that a™ term can be neglected, then the above 
expressions provide the particle's orbit elements in a drag-dominated ring, which are 

e(x) ^ — ^^^^ and smmuj{x) ^ '^^ (98) 

where (i(x) ^ x (see Eqn. [HS]) and with q = C^/ (m — e). The satellite's radial torque density 
is Eqn. ( I85l) . so the total torque that the satellite exerts on the ring is the integral 

Ts = —^da = mnaooa^'^^ / essm{mLu)dx' (99a) 

J ring fJUi J j-jng 

= -emTraooa^^T^sQ / ^ = ?1| = ^gt (99b) 

when it is assumed that the ring has a constant surface density {s = 1) that extends 
everywhere. But if the ring instead had a sharp edge at the resonance, then the satellite's 
torque would be half that. 



The integrand in Eqn. (]99bp indicates that the satellite's torque on the ring is exerted 



over a fractional radial scale |x| torque ~ which makes sense since this is where the stream- 
lines' e and \u\ are maximal (Eqns. 1^51) . But keep in mind that Eqns. (1981) also ignored 
the ring's internal forces a^, which Section 12.4.21 showed to be important in the nonlinear 
zone whose radial extent is a; at/, = a/I^'sT from the resonance (Eqn. HI]). Those internal ring 
forces also tend to inhibit streamline-crossing, which they achieve by reducing the stream- 
line's eccentricities below that given in Eqn. fl98l) . This in turn reduces the satellite torque 
to something less than Eqn. ( I99b[) . But this torque reduction should still be insignificant 



whenever xnl ^ |a;|torque, when the revised drag parameter satisfies Cd ^ 

however, this criterion is not satisfied, then the Goldreich-Tremaine torque formula, Eqn. 

(!95|) . will overestimate the satellite's torque on the ring. 

Of course, the dissipation in the ring models of Sections I3.2.1H^TT^ and Figs. [TH5] was 
due to viscosity rather than drag. However, the preceding discussion suggests that Eqn. (l95l) 
will provide a reliable estimate of the satellite's torque whenever Aa^, which is defined here 
as the radial distance over which the satellite substantially excites the streamlines' uj, satisfies 
Aa^ ^ Aa^L where AajvL = \/\ips\0'r is the width of the nonlinear zone in physical units. 
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Note that the width of Mimas' m = 2 ILR is Aa^L ^ 20 km (Section [232]), while Figs. HHS] 
show that these viscous ring models have Aa^ < 10 km. Consequently, the satellite's torque 
is reduced below Eqn. f P^ . which explains those models' difficulty in achieving a torque 
balance at the ring's outer edge. 

The torque that the satellite exerts on each model ring is also reported in all of the 
figure captions. Those captions show that the ring's internal forces — gravity and pressure — 
reduce the satellite's torque below Eqn. (^5^ by three to seven orders of magnitude. This 
explains why the viscous ring model of Fig. [7] requires a viscosity ratio z/f,/z/s ~ 10^, which 
also boosts the satellite's torque by that factor, in order to achieve a torque-balance at the 
B ring's outer edge. The exception of course is the drag-dominated ring of Fig. [91 In that 
model, the satellite's torque is only about three times smaller than Tqt, with the greater 
torque efficiency being due to the drag's ability to excite larger u across a wider radial span 
in the ring such that Aa^ > Aqnl- 



3.2.8. assessing the ring's internal forces 

The relative importance of the various accelerations that a ring particle experiences due 
to ring gravity, pressure, viscosity, etc., is assessed by dividing the real part of the complex 
equation of motion f[5T]) by the satellite's forcing function t/^s, which yields 

{Acp + -Ag + Ap) cos(mcj) - {Ay + Ad) sin(mu;) + 1 ~ (100) 

where Acp = ede/ips, Ag = ^e{a'^^)/ips, Ap = ^e{a'^^)/i/js, Ay = ^e{ipu)/i/js, and Ad = 
5^e('(/'drag)/'^s, where ipu and V^drag are Eqn. and (IMl) . Note that Eqn. fllOOp is only 
approximately true because the accelerations appearing in the equation of motion (ISTi) are 
complex, yet Eqn. fllOOl) neglects their smaller imaginary parts. The quantities Ag, Ap, 
Av, and Ad represent the strength of the m^^ component of a particle's acceleration due 
to ring gravity, pressure, viscosity, and drag, all in units of the satellite's forcing ifjs, while 
Acp = eed/ips is the relative strength of the centrifugal and Coriolis accelerations that the 
central planet exerts on the particle due to its noncircular motions. 

The relative accelerations A are shown in Fig. [11] for the two models that achieved a 
torque-balance at the ring's outer edge: the model described in Section 13.2.51 and Figs. U\M 
that invokes an extreme viscosity ratio i^b/i^s ^ 1 to achieve its torque-balance (see Fig. 
[TTh). and the model of Section [3.2.61 and Figs. [^HTUlthat relies on drag to balance the ring- 
satellite torques (Fig. [TTb). Both models show that the B ring's internal accelerations are 
small interior to the ring's outermost ~ 40 km. There, the particle's motion balances the 
satellite's forcing against the central planet's centrifugal/ Coriolis forces, i.e., Acp — —1 since 
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|u)| <^ 1, which is equivalent to the single particle solution, Eqn. fl29al) . Figure [TT] also shows 
that, in the model B ring's outermost ~ 40 km, the ring's gravity is the dominant internal 
ring-force whose strength is comparable to the satellite's forcing and the central planet's 
centrifugal/ Coriolis forces. Those curves show that the acceleration due to ring pressure is 
small everywhere except at the ring's outer edge, whose effects there are only marginally 
resolved in these models (see Section 13.2.31) . Figure [11] also shows that viscosity has no 
direct effect on the ring's epicyclic amplitude; instead, its influence enters indirectly via the 
torque-balance Eqn. (1551) . Figure [TTb also illustrates how strong the drag force must be if it 
is indeed responsible for a torque-balance at the B ring's outer edge. 

Lastly, recall Eqn. fl82|) . which shows that the m = component of the ring's conser- 
vative accelerations, + a^^, can displace the location of the resonance, which is the site 
where d{x) = 0. It turns out that this displacement is quite small — too small to be resolved 
in the models considered here. However, that displacement is easily inferred from a linear 
interpolation of the d{xi) data that the streamline model generates. For instance, an inter- 
polation of the data generated by the two B ring models in Figure [11] shows that the ring's 
internal forces displaces Mimas' m = 2 ILR inwards about 15m. That small displacement is 
due to ring gravity, since the ring pressure is negligible there. 



4. Discussion 



4.1. Heating the ring's outer edge 



Another issue that me rits consideration is the viscous heating of the ring's outer edge 



(e.g, iBorderies et al.lll982l ). which might be important when collisions are the dominant 
source of viscosity, since t/ oc in this case. To assess this heating, one has to calculate 
the rate at which the ring's dissipative forces do work on each ring particle. That quantity 
is probably positive, because the viscous delivery of orbital energy to the ring's outer edge 
likely exceeds the rate at which the satellite with draws orbital energy f rom the ring edge, 
resulting in a dynamical heating of the ring-edge (IBorderies et al.lll982l ). Of course, other 
processes also tend to cool the ring particles' random velocities, such as dissipative collisions 
that convert impact energies into thermal heat, mechanical grinding and fracturing of ring 
particles, and thermal radiation. Those other cooling mechanisms have not yet been consid- 
ered, but do need to be included in the ring's energy-balance equation, which could be used 
to relate the ring's viscosity z/(a) to the ring particles' dispersion velocity c{a). But many 
of the terms in that equation will be difficult to quantify, and will be deferred to a followup 
study. One might have to resort to an order-of-magnitude type analysis, or perhaps settle 
for an upper limit on the heating that occurs at the ring's edge. However, upcoming Cassini 



-44- 



observations of the rings during the August 2009 equinox will be helpful here, since that is 
when the Sun will pass through the ring plane. If the ring-edge is dynamically hot and thick, 
then the shadows it will cast onto the ring-plane should be observable, which would allow 
one to infer the ring-edge's vertical scale height h and the ring particles' dispersion velocity 
c = hVL there. 



4.2. Future applications 



The outer edge of Saturn's A ring is maintaine d by an m = 7 ILR with the coorbital 
satellites Jan us and Epimetheus (IPorco et al.lll984j ). whose semimajor axes differ only by 
about 0.03% ( Jacobson et al. 2008h. These satelli te's mutual attractions cause their orbits 
to swap about every four years (lYoder et al.lll983l ). such th at only one coorbit al appears to 
have its m = 7 ILR in the A ring at any instant of time (jSpitale et al.ll2009l ). One might 
be tempted to use this streamline model to calculate the ring's response to each individual 
satellite, and then to superimpose the model outputs to obtain the ring's total response 
to both satellites. However this might not lead to reliable results, due to the satellites' 
time-varying orbits. Because the torque that they exert on the A ring changes periodically 
with time, the outer A ring never experiences a balance of the viscous and satellite torques 
at any given instant, and so the model's torque balance Eqn. fl5^ does not apply here. 
Instead, a time-average of the outer A ring's viscous torque should be balanced against 
the time-averaged torque that the coorbitals exert on the ring. However, the calculation of 
those time-averaged torques is subtle, and will be saved for a foUowup study of the A ring 
JSpitale et al.ll2009h . 



The streamline model can also be used to study narrow eccentric ringlets, but with some 
revision. For instance, if ringlets are maintained by small unseen shepherd satellites, then 
the impulse approximation should be used to calculate the torque that they exert on the 
nearby ringlet. Also, the differential precession that is due to planetary oblateness, which is 
important in a ringlet system, will also need to be accounted for in a revised version of the 
streamline model. 

Lastly, the streamline model can also be modified so that it can simulate linear as well 
as nonlinear spiral density waves. Spiral waves transport angular momentum through a ring 
via the tangential accelerations that particles experience due to ring gravity and pressure. 
Those terms are negligible in this study of a nearly peri-aligned B ring, but their inclusion 
will be needed in order to handle the spiral waves' angular momentum transport. Those and 
other related problems will be considered in followup studies of planetary rings. 
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5. Summary of Results 

The preceding describes a model of a broad, sharp-edged planetary ring that is confined 
by a satellite's r ri^^ L i ndbla d resonance. This model utilizes the streamline formalism of 



Borderies et al.l (Il982l . Il985l ). which makes the calculation of the ring's internal forces — 
ring gravity, pressure, and viscosity — quite tractable. The model also includes a simple 
prescription for handling the drag forces, such as such as plasma, Yarkovsky, atmospheric, 
and/or PR drag, that small ring particles might experience. The model's main inputs are 
the ring's surface density doo, the ring particle's dispersion velocity c, the ring's kinematic 
shear and bulk viscosities z/^ and u^, and a dimensionless drag coefficient Cd- This streamline 
model solves a nonlinear form of Newton's second law of motion to obtain the streamlines' 
orbital eccentricities e(a) and longitudes of perihelia a) (a) as functions of the streamlines' 
semimajor axis a. The model also balances the ring's viscous torque against the satellite's 
gravitational torque in order to calculate the ring's surface density a{r, 9) as a function of 
radius r and longitude 6. That analysis also shows how to use linear and angular momentum 
ffuxes to calculate the effects of viscosity and pressure, both of which are discontinuous at a 
ring's sharp edge. 

The streamline model is then applied to the outer edge of Saturn's B ring, which is 
maintained by an m = 2 ILR with the satellite Mimas. A suite of B ring scenarios are 
examined in order to illustrate how the model outcomes depend upon the the ring's physical 
properties (Too, c, z/^ and Vb, and possibly C^, with the main ffndings listed below. 

1. As one might expect, increasing the ring's surface density tends to decrease the 
ring-edge's epicyclic amplitude Rm- 

2. Pressure in the ring is controlled by the ring particle's dispersion velocity c, but 
increases in c are manifest only at the ring's outer edge, which is where the ring's pressure 
drop is greatest. However the consequences of that pressure drop are not fully resolved in 
the models considered here, which have radial samplings of 0.5 to 2 km. 

3. The magnitude of the B ring's internal forces are compared, and it is shown that the 
ring's gravity dominates over its other internal forces that are due to pressure and viscosity. 
In the B ring's outer ~ 40 km, the gravitational force that the ring exerts on a particle is 
comparable to the satellite's forcing. However, at sites well interior to the ring's outer ~ 40 
km, ring gravity is relatively small, so the ring particles adopt the familiar single-particle 
solution e = \iIJs/x\ there. 

4. The ring's viscous torque is controlled by its shear viscosity Vs-, which also governs 
the rate of the ring's radial spreading. However the satellite's torque on a viscous ring, which 
opposes that spreading, is proportional to the angle by which the streamline's periapse lags 
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behind the satelhte's longitude. That lag angle is a linear combination of the ring's shear 
and bulk ui, viscosities, so the satellite's torque on the ring is sensitive to the sum of 
those viscosities. Interestingly, a conventional ring model that has a bulk viscosity that is 
comparable or less than the shear viscosity fails, by a very wide margin, to balance the 
satellite's torque against the ring's viscous torque, so that scenario fails to account for the B 
ring's sharp edge near a resonance. Nonetheless, increasing z/f, does strengthen the satellite's 
torque on the ring, and models show that a torque balance becomes possible if I'b/^s ~ 10"^. 
In other words, the ring particles' shearing tangential motions must be relatively free of 
frictional dissipation, while the ring particles' compressed radial motions must result in heavy 
frictional dissipation. Such a scenario might be possible if particles in the B ring's compressed 
regions are confined shoulder-to-shoulder so that there is little voidspace between the ring 
particles there, and ring's volumetric density becomes incompressible there. In this case, 
ring particles must rise and then fall vertically, sliding or perhaps tumbling past one another 
as the satellite's perturbations drives the ring-edge radially inwards and then outwards. This 
periodic avalanche of ring particles could be quite lossy, possibly resulting in a very large 
bulk viscosity Uh. However it is unclear whether the B ring edge actually behaves in this 
way, and whether its viscosity can satisfy this remarkable requirement of I'b/^s ~ 10^, so this 
finding is speculative. Note also that current ring observations only provide measurements of 
the combined viscosity Ug + fvb where / ~ 0{1), so the ratio Vb/vs in a perturbed planetary 
ring-edge is actually unknown. 

5. Drag forces provide an alternate means of boosting the satellite's torque on the ring. 
A generic drag force is considered, one that is proportional to a ring particle's noncircular 
velocity, and its main effect is to damp the particle's eccentricity while causing its longitude 
of perihelia u) to trail behind the satellite's longitude. Since the satellite's torque on the 
particle is proportional to sin(ma)), this allows the satellite to torque a wide annulus in the 
ring, which also makes a torque-balance at the ring-edge favorable. However this torque- 
balance due to drag can only be effective if the ring particles' mass distribution is dominated 
by small particles, which is probably not the case for Saturn's B ring. 

6. The torque that the satellite exerts on the entire ring is calculated for a variety of 
viscous B ring models, and it is shown to be 3-7 orders of r nagnitude smaller than that 



anticipated by the Goldreich-Tremaine torque formula Tqt (iGoldreich &: Tremaind Il978 



19821 ). This is due to the ring's self-gravity, which suppresses the streamline's eccentricities 
in the nonlinear zone that is near the resonance, which also reduces the satellite's torque 
there. However, a drag-dominated ring can experience a much larger torque. This is due to 
the streamlines' trailing longitudes of perihelia a), which communicates the satellite's torque 
beyond the nonlinear zone, and can result in a torque that is comparable to Tqt- 
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7. The outer edge of the B ring could he as far as Aa ~ 24 km beyond Mimas' m = 2 
inner Lindblad resonance, which is the site where the ring's viscous torque precisely balances 
the satellite's torque on the ring. Models of the B ring are adjusted so that the simulated 
ring-edge's epicyclic amplitude agrees with the observed Rm = 45 km amplitude, and that 
a torque balance is achieved at distances of < Aa < 24 km beyond Mimas m = 2 ILR. 
This requires the outer B ring to have a surface density of 10 < doo < 280 gm/cm^ in 
the ring's outermost ~ 40 km, with the larger surface densities required for rings whose 
edges lie farther beyond the resonance. A more detailed comparison of models to Cassini 
observations, which is ongoing, will lead to a more precise measurement of that ring's surface 
density. Note, though, that the displacement of the ring-edge from the resonance is not due 
to the ring's internal forces altering the resonance location. For instance, the displacement 
of the resonance due to ring self-gravity and pressure is tiny, only about 15m for the models 
considered here. But if Mimas had instead been more or less massive, then the B ring's outer 
edge would then lie interior or exterior to its present location. 

8. Models of the B ring predict that the ring-edge's surface density should be enhanced 
by ~ 50% at the ring's longitude of peripase. The ring's outer edge should also exhibit a low 
surface density shoulder at the ring's longitude of apoapse. These surface density variations 
are due to the satellite's perturbations, which compresses the streamlines at periapse and 
rarefies them at apoapse. It is curious that these periapse enhancements have not been 
reported in spacecraft observations of the B ring edge, but this non-detection may be due to 
a saturation of the ring's surface brightness I/F t hat is expected to occur when the ring's 



optical depth exceeds about 0.3 (jPorco et al.ll2008l ) 



9. Modifications to the streamline model will also allow its application to other dense 
planetary rings, such as the outer edge of Saturn's A ring, and the many narrow and some- 
times eccentric ringlets that orbit both Saturn and Uranus. And with additional physics, 
this model will also provide a useful tool that can be used to simulate nonlinear spiral den- 
sity waves. Detailed comparisoii s of models to spacecraft observations of Saturn's rings are 



ongoing (c./. ISpitale et al.ll2009l ). and that activity should yield better estimates of, or else 
place limits on, the ring's physical parameters (Too, c, and Uh, and Cd- Such studies will 
also lead to a better understanding of the mutual interactions that are exerted in these very 
interesting ring-satellite systems. 
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Table 1. Physical properties of the Saturn-Mimas system'^ 



Saturn's GMp 


3.79312077 x 10^ km^/sec^ 


Saturn's zonal harmonic J2 


1.629071 X 10-2 


Ja 


-9.3583 X 10-^ 


Je 


8.614 X 10-5 


Js 


-1.0 X 10-5 


Mimas' angular velocity^ Qs 


381.9944522 degree/day 


Mimas' fractional mass Hs 


6.5969x10-^ 



'^from iJacobson et al.l (120061 ) , except where noted otherwise 
^as reported by the JPL Solar System Dynamics website 



http://ssd.jpl.nasa.gov/7sat_elem on June 26, 2008 
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Fig. 1. — Simulations of the B ring's epicyclic amplitude Rm are plotted versus a — ar, which 
is the radial distance from Mimas' nominal m = 2 ILR. Solid curves are for simulated rings 
that have a variety of surface densities cxoo that are indicated above in units of gm/cm^. The 
dashed cXoo = curve is the single particle solution, Eqn. (l28ll . and the h orizontal dotted 
line i ndicates the B ring's observed epicyclic amplitude of i?obs — 45 km (ISpitale fc Porco 



20061 ) ■ The model rings all have viscosities of Vg = i^b = 50 cm^/sec, a stability parameter of 



Q = 2, and no drag (C^ = 0.) The rings' outer edges are also forced to reside at the nominal 
resonance where x = = a — a^; see Section 13.2.11 for details. All of these simulations 
use = 300 streamlines that are distributed uniformly at the B ring's outer edge. The 
streamlines' radial widths differ in each model, with Aa = 0.50,0.50,0.67, 1.0, and 2.0 km, 
with wider streamlines used in the heavier rings for reasons explained in Section 13.11 The 
total radial extent of the modeled region is w = NAa = 150 km in the lightest ring and 
w = 600 km in the heaviest ring. The torque Tg that the satellite exerts on the ring is Eqn. 
(!86|) evaluated at the ring's outer edge, and it is quoted in units of Tgt (Eqn. [95]). Those 
torque ratios range over \Ts/Tgt\ = 7.6 x 10~^ for the lighter (Too = 10 gm/cm^ ring to 
\Ts/Tgt\ = 1-7 X 10~^ in the heavier a^o = 1000 gm/cm^ ring. 
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Fig. 2. — The nonlinear parameter q is plotted versus distance from resonance for the B 
ring models of Fig. [H 
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Fig. 3. — Results of three simulations of the outer B ring as it is perturbed by Mimas' m = 2 
ILR. The three simulated rings have outer edges that lie at radial distances of a — a,. = 0, 10, 
and 25 km beyond the nominal m = 2 resonance position. Labels indicate each ring's surface 
density in units of gm/cm^, which is adjusted so that the epicyclic amplitude Rm at the 
ring's outer edge is comparable to the observed value -Robs — 45 km, which is indicated by 
the dotted line. These models use = 300 streamlines having widths of Aa = 0.50, 0.87, 
and 1.42 km, with wider streamlines being used in the heavier rings, so the total radial 
extent of the modeled regions are 150, 260, and 425 km. The remaining model parameters 
are identical to that used in Figs. [THll with I's = = 50 cm^/sec, Q = 2, and Cd = 0. The 
torque that the satellite exerts on the ring ranges over \Ts/Tgt\ = 7.6 X 10~^ in the lighter 
(Too = 10 gm/cm^ ring to \Ts/Tgt\ = 2.9 x 10~^ in the heavier cXoo = 280 gm/cm^ ring. 
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Fig. 4. — The B ring's epicyclic amplitude Rjn{a), its fractional surface density s{a), and 
its longitude of periapse a;(a) are plotted versus distance from resonance a — ar for models 
that all have the same surface density (Too = 30 gm/cm^, viscosities Us — ui, — 50 cm^/sec, 

and stability parameters of Q = (which corresponds to a pressureless ring, black curve), 
Q = 1 (yellow), Q = 5 (cyan), Q = 10 (green), Q = 11 (red), and Q = 12 (blue). The Q = 
to 5 simulations are nearly indistinguishable and lie under the black curves. There is no 
drag in these simulations (C,i = 0), and the rings' outer edges are at the nominal resonance 
where a = a^.. These simulations use = 300 streamlines that are distributed uniformly 
over the ring's outermost 150 km, so the radial resolution here is 0.5 km. The torque that 
the satellite exerts on the ring ranges is \Ts/Tgt\ — 1-5 x 10~^ in the < Q < 5 models, 
7.7 X 10-5 in the Q = 10 model, 6.1 x 10"^ for Q = 11, and 2.3 x 10"^ for Q = 12. 
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Fig. 5. — The simulations of the B ring's epicychc amphtude Rm (Fig. a), fractional surface 
density s (Fig. b), and longitude of periapse u (Fig. c) are all plotted against the radial 
distance from resonance a — Ur- Figure d plots the rings' dimensionless angular momentum 
luminosity due to viscosity £1, and the satellite's gravitational torque is, as well as the rings 
nonlinearity parameter q. Three simulations are shown, and they all have an undisturbed 
surface density o"oo = 30 gm/cm^, stability parameter Q = 2 (which corresponds to particle 
dispersion velocity of c = 0.82 mm/sec and a ring vertical half-thickncss oi h = 5.4 meters). 
The ring's outer edge is also placed at resonance at a = a^, and the drag coefficient Cd = 0. 

= 300 streamlines are used to model the ring's outermost w = NAa = 150 km with 
a spatial sampling of Aa = 0.5 km. These three simulations do have distinct viscosities 
I's — i^b — 5, 50, and 500 cm^/sec, which are indicated by the labels in Fig. c. Note also that 
the curves in Figs, a, b, and d all lie on top of each other for each of the three simulations, 
due to Rm, s, q, and i being insensitive to the choice of i/g — i^b- The torque that the satellite 
exerts on the ring is \Ts/Tgt\ = 1-5 x 10~^, 1.5 x 10~^, and 1.5 x 10~^ for the i^s — J^b — 5, 50, 
and 500 cm^/sec ring models. 
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Fig. 6. — The cpicyclic amplitude Rm (Fig. a) is plotted versus radial distance a — from 
the simulated B ring's outer edge, as well as its fractional surface density s (Fig. b), longitude 
of periapse Co (Fig. c), and the angular momentum luminosities 1^ and Is (Fig. d). Shown are 
the results of three simulations that all have the same undisturbed surface density cxoo = 30 
gm/cm^, stability parameter Q = 2, no drag (C^ = 0), and an edge at the nominal resonance 
where a — ar- N — 300 streamlines are used to model the ring's outermost w — 150 km with 
a spatial samphng of Aa = 0.5 km. The simulated rings' shear and bulk viscosities all satisfy 
i^s + T-'b = 100 cm^/sec while having distinct ratios fh/^s = (green curve), 1 (red), 1000 
(blue), and 7873 (black curve, which is the only simulation that satisfies the torque-balance 
requirement). Note that simulations having Ub/'^s = and Uh/us = 1 are indistinguishable in 
Figs. a,b, and d, so the green curves are hidden under the red. The torque that the satellite 
exerts on these model rings ranges over 8.8 x 10~^ < \Ts/Tgt\ = 1-5 x 10^^. 
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Fig. 7. — Epicyclic amplitude Rm (Fig. a), fractional surface density s and nonlinearity pa- 
rameter q (Fig. b), longitude of periapse a) (Fig. c), and the angular momentum luminosities 
£1, and £s (Fig. d) are all plotted versus radial distance a — ar for a simulated B ring that 
has an undisturbed surface density (Too = 226 gm/cm^, stability parameter Q = 2, a kine- 
matic shear viscosity of Ug = 0.00603 cm^/sec, kinematic bulk viscosity = 100 cm^/sec, 
no drag {Cd = 0), and with an outer edge that lies 24 km beyond the nominal resonance 
position. N = 300 streamlines are used over the ring's outer 424 km, so this calculation has 
a spatial sampling of Aa = 1.41 km. The torque that the satellite exerts on the model ring 
is \Ts/Tgt\ = 4.2 X 10-^ 
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Fig. 8. — Equation fl93ap is used to calculate the ring's relative surface density ojoao as a 
function of semimajor axis a for the viscous B ring model that is described in Fig. [71 Equation 
(I93bp is then used to convert that surface density profile into a function of planetocentric 
distance r, which is plotted here versus distance Ar = r — from the nominal resonance. 
These curves give the simulated B ring's relative surface density along longitude 9 = 9s 
(which is along the ring's longitude of periapse), longitude 6* = 6*^ ± 90° (along the ring's 
longitude of apoapse), and along the intermediate longitude 9 = 9s^ 45°. 
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Fig. 9. — This simulated B ring has a surface density — 265 gm/cm^, stabihty parameter 
Q = 2, a kinematic shear viscosity of = 46 cm^/sec, kinematic bulk viscosity = 0, and 
a drag parameter = 1.0 x 10^"', with these parameters chosen so that epicyclic amplitude 
at the ring's outer edge is i?„t = 45 km, and that the ring's viscous torque balances the 
satellite's torque at a distance of a — = 24 km beyond the ring's nominal resonance. 
N = 300 streamlines were used to model the ring's outermost 424 km, so the spatial sampling 
here is Aa = 1.41km. The torque that the satellite exerts on the ring is \Ts/Tgt\ = 0.32. 
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Fig. 10. — Equations ( l93i) are used to calculate the ring's relative surface density cr(r, 0)/aoo 
as a function of radial distance Ar = r — from the nominal resonance for the model B ring 
that is described in Fig. [91 These curves give the simulated ring's relative surface density 
along the ring-edge's longitude of periapse, its longitude of apoapse, and at an intermediate 
longitude. Note also the small bumps seen at the rightmost part of these curves. They are 
due to the small, marginally resolved surface density excess that is barely seen to the right 
in Fig. [9b. That bump is due to the pressure drop that the outermost streamline experiences 
(e.g, Eqn. [57]), and it disappears when the ring is pressureless with Q = 0. 
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Fig. 11. — These figures plot the relative accelerations that a ring particle experiences due 
to ring gravity Ag (black curve), pressure Ap (green curve), viscosity Au (red curve), drag 
Ad (orange curve), and Acp, which is the acceleration that the central planet exerts on the 
particle; see Section 13.2.81 for details. The upper figure is for the model reported in Figs. 
[THHlthat achieved its torque-balance at the ring's outer edge via an extreme viscosity ratio 
i^b/'^s = 8473. The lower Figure is for the model described in Figs. [5HTUI which relies on drag 
to enable a torque-balance. 



